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THE RECEPTIVITY OF BOUNDARY LAYERS ON BLUNT BODIES 
TO OSCILLATIONS IN THE FREE STREAM 

By 

Chester E, Grosch* 

INTRODUCTION 

The prediction of the position of transition on s body is one of the 
most difficult problems in fluid dynamics. Linear stability theory with the 
H e n " criterion (see Mack, 1977, for an extended discussion) is probably 
the most widely used method of predicting transition. The major weakness of 
this method is that it employes a relative measure of flow instability; the 
growth of some measure of the disturbance size, whether amplitude, energy, 
or whatever, relative to an unknown initial size of the disturbance. 

It is clear, as pointed out by Mack (1977), that the position of tran- 
sition must depend on the absolute size of the disturbance, and thus on the 
initial size of the disturbance as well as the growth rate of the unstable 
disturbance. If there were no disturbance at all, that is the initial size 

were zero, the stability or instability of the flow would be irrelevant. 

* 

The flow would remain laminar, and transition would not occur anywhere on 
the body unless, of course, the flow separated, On the other hand, if the 
initial disturbance were very large, say of the order of the mean flow 
speed, then transition would, presumably, occur very near the front of the 
body, perhaps in the immediate vicinity of the stagnation point or line. 

A rough estimate of the magnitude of the initial amplitude of the 
TollmienrSchlichting waves at the beginning of the region of instability can 
be made if it is assumed that the M e n " criterion is an approximately valid 
transition method and if one accepts the estimate of Klebanoff, Tidstron, 
and Sargent (1962) that, at the beginning of transition, the rms value of 
the perturbation velocity is about 1.5% of the free stream speed. Since the 


♦Professor, Department of Oceanography, Old Dominion University, Norfolk, 
Virginia 23508, 


initial proposal of th« "« n " criterion by Smith (1956) and Van Ingen 
(1956), extensive studies of transition data hava shown that n varies from 
about 2,5 to 10.0, depending on the test facility. Taking the disturbance 
amplitude at transition to be 0.015 of the free stream speed, and n to 
vary from 2.5 to 10,0, it follows that the amplitude of the disturbance at 
the beginning of the region of flow instability is in the range 10 -3 to i0 -6 
of the free stream speed. This suggests that a very weak coupling of free 
stream disturbances to Tollmien-Schlichting waves in the boundary layer may 
be extremely important in the initiation of th* transition process. If the 
free stream disturbances have amplitudes of O(c), even disturbances in the 
boundary layer of 0(e 2 ) may play a role in initiating transition. 

In order to incorporate the influence of the flow environment in sta- 
bility prediction methods, and thereby provide a rational basis for theoe 
methods, it is necessary to be cble to calculate the initial size of the 
disburbance caused by external, i.e, free stream, disturbances. There are 
four classes of free stream disturbances which may be responsible for the 
generation of the Tollmien-Schlichting waves which lead to transition 
(Obremski, Morkovin, and Landhal, 1969). They are: mean flow unsteadiness, 

free stream vor'cicity, sound, and entropy fluctuations, If the fluid were 
truly incompressible, only the unsteadiness of the mean flow and free stream 
vorticity could exist and initiate the transition process, In reality, 
sound waves are always present. Entropy fluctuations appear to occur only 
at high supersonic speeds. In a recent study, Harvey and Bobbitt (1981) 
have examined the anomolies between wind tunnel and flight transition mea- 
surements. They state that the experimental results show"... the dominance 
of free- stream disturbance level on the transition process from beginning to 
end." Further, they conclude that there is a question of"... whether sound, 
unsteadiness, or spectral peaks is the most influential on boundary layer 
receptivity,..." Morkovin (1978) has argued that, in medium to low subsonic 
flows, sound waves can be modeled quite accurately by a time varying 
oscillation superimposed upon a steady mean flow, provided that the 
wavelength of the sound wave is much larger than the characteristic length 
scale of the body. For a sound wave impinging on the nose of a body this 
characteristic scale length would appear to be the radius of curvature of 
the body at the forward stagnation point. 
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Recent studies, (Salwen, Kelly, end Grosch, 1980; Grosch end Salwen, 
1980; Selwen and Grosch, 1981) heve been addressed to the problem of calcu- 
lating the amplitudes of the discrete end continuum eigenfunctions of the 
linearized, parallel flow stability problem. The results of these studies 
can be sunmarized as follows: first, the proper inner products for ortho- 

gonality of these eigenfunctions were found; second, it »“'j shown that the 
set of discrete and continuum eigenfunctions for both the temporal and spa- 
tial problems is complete; and third, a procedure for calculating the ampli- 
tudes of the discrete and continuum eigenfunctions was found. In order to 
calculate these amplitudes, the disturbance must be known, as a function of 
time, on a plane perpendicular to the boundary. 

In order for this procedure to be applied it is necessary for the dis- 
turbance to be small so that linearization is valid; that the parallel flow 
approximation be valid, at least approximately; and that the form of the 
disturbance be known throughout the boundary layer and free stream at one 
location on a two-dimensional body. Even if the disturbances are small, 
this theory cannot be applied in the immediate vicinity of the leading edge 
of a body because the flow in that region is not even approximately 
parallel. A basic problem is to determine the flow within the boundary 
layer caused by disturbances propagating and being advected by the mean flow 
toward the front of a blunt body, 

The response of the boundary layer on a body to oscillations in the 
free stream has been the subject of a number of studies beginning in the 
early 1950's. Some representative examples are the work of Moore (1951), 
Lighthill (1954), Stuart (1955), Rott (1956), Lin (1956), Glauert (1956), 
Carrier and Di Prima (1957), Gibson (1957), Watson (1959), Lam and Rott 
(1960), Rott and Rosanzweig (1960), Sarma (1964), Stuart (1966), and 
Ackerberg and Phillips (1972). Riley (1975) and Schlichting (1979, Chapter 
15) give comprehensive reviews of recent work in unsteady boundary layer 
theory. None of these authors have considered the application of their 
results to the boundary layer receptivity problem. 

These studies have considered either the unsteady flow in the vicinity 
of a stagnation point or on a flat plate. The flat plate problem presents 



great difficulties because of the singularity of the boundary layer equa- 
tions at the leading edge of the plate. Carrier and Di Prima (1957) studied 
this problem within the framework of a linearized, Oseen approximation, to 
first order in the amplitude of the oscillation, ‘they found that, far from 
the leading edge of the plate, the time dependent part of the flow is a 
Stokes shear wave. Ackerberg and Phillips (1972) studied the same problem 
using asymptotic, again to first order in the amplitude of the oscillation, 
and numerical techniques. They also found that, far downstream of the lead- 
ing edge, the flow develops into a nearly Blasius steady flow and a Stokes 
shear layer imbedded within the steady boundary layer. Htte most surprising, 
and unexpected, result of Ackerberg and Phillips is that "... most flow 
quantities approach their asymptotic values far downstream through damped 
oscillations." 

Quite recently, Grosch and Salwen (1982) have found a solution of the 
Navier-Stokes equations for an incompressible stagnation point flow whose 
magnitude oscillates in time about a constant, nonzero, value (see Attach- 
ment 1). The analytic solutions, to third order in the amplitude, obtained 
using inner and outer expansions, are in complete agreement with the results 
of numerical integrations. The mean flow is the steady stagnation point 
flow plus second, and higher, order flows driven by the Reynolds stress of 
the oscillatory components. The oscillatory flow is that of the fundamental 
and all of the higher harmonics of the fundamental. 

The mean flow is, to lowest order, independent of the oscillating flow. 
The structure of the oscillatory components of the flow depends on the fre- 
quency of oscillation, not on the distance from the stagnation point. As 
was expected, the oscillatory flow develops into a Stokes shear layer at 
high frequencies. The results of Carrier and Di Prima and Ackerberg and 
Phillips that the Stokes shear layer develops far downstream may be due to 
the use of the boundary layer equations in a region in which they are 
invalid. 

In any case, the flat plate with zero thickness and a "sharp" leading 
edge is a nonphysical model; all physical bodies have a blunt leading edge, 
at least if a continuum approximation can be applied to the flow past the 
body . 
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In this paper we examine one aspect of the boundary layer receptivity 
problem; that of the flow induced in the boundary layer on the forward por- 
tion of a blunt body by unsteadiness in the mean flow. As was mentioned 
above, this is also an approximate solution for the disturbances generated 
in the boundary layer by a long wavelength sound wave impinging normally on 
the front of a body in steady motion. This work is an extension of our 
previous study of the flow in the neighborhood of a stagnation point (Grosch 
and Salwen, 1982). 


FORMULATION 

Consider the flow past a blunt, two-dimensional body. The usual bound- 

fsi 

ary layer coordinates ate used! x is the distance along the surface of the 
body from the forward stagnation point and y is the distance from the body 

IV 

along the normal to the body. A solution is to be found for x 20 of the 
time-dependent, two-dimensional boundary layer aquations 


aw FV* 

3 x 3y 


3 u 

3 u 

3u 

3U 


T u + V — , * 

r* a# 

, „ 

+ U — 

- A* 

3 t 

3x 

3y 

3 t 

3x 



( 1 ) 

( 2 ) 


with the boundary conditions 

u(x,o,t) * v(x,o,t) - 0, (3) 

u(x,y,t) + U(x,t) as y + ». (A) 

A# 

Here, (u,v) are the components of the velocity in the x and y 
directions, U(x,t) is the component of the potential flow which is tan- 
gential to the body surface, and v is the kinematic viscosity. It will be 
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ft » 


c 

l 

I 

assumed chat the potential flow is a temporally oscillating flow, super- 


imposed on a steady mean flow; specifically it is assumed that, l 

't 

U(x,t) « U H<x/£> [l + e cos tot], (5) i 

W I' 

where U 0 is Che magnitude of the mean potential flow, | 

£ 

l is an appropriate length scale of the body, ? 

t is the dimensionless amplitude of the oscillating component, ■ 

r 

and , i 


H is a dimensionless function, giving the variation along the body 
of the tangential component of the potential flow. 

It will be assumed that the body is symmetric and is at a zero angle of 
attack, so that H is of the form 

n<*/t> • ? , 2ntI <;/«> 2 " +1 , (6) 

n"0 

with the {aj} dimensionless constants. Specific examples will be 
considered below. 

METHOD OF SOLUTION 
Define dimensionless variables 


K - x/A, 

h * (ajRe) l ^ 2 (y/l) 

T - aiU Q t/A, 

and parameters, a Reynolds number 

Re - U t / v 
o 


(7) 

( 8 ) 
(9) 


( 10 ) 


| 

f 
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and a Strouhal number, a dimenaionlesa frequency, 

o * uit/(ajU o ) , (u) 

Defining a stream function auch that 

u w /3 y , v ■ -3i|»/3x, 

by 

♦ " [<U o £)/(a v Re) l/2 ] ♦(£ ,H ,t), 
equation (1) ia satisfied identically and (2) is 

3 3 4> - l fa$ 3 2 $ 3* 3 2 * 1 3 2 $ „ /tx . 

— 2L + ai I-— — - -I I— I— « eoH(5) sin or 

3n 3 [H 3n2 3n 353ti J 3t3n 

- a* H(0H'(0[ i + e cos ot] 2 , 


(12a, b) 
(13) 


The boundary conditions, (3) and (4) are 


34 34 

JL. . JL « 0 

3n * 35 


at n ■ 0, 


and 


(15) 


+ H(5)[l + e cos or] a$ n ♦ (16) 

3n 

The tangential component of potential velocity, H(5), is a power series in 
?. It is clear, from (14), that $ must also be a power aeries in 
with each coefficient a function of h, x, and the parameter e. Each of 
these coefficients wilt be expanded in a Fourier series in t, with each 
Fourier amplitude a function of n and the parameter e. Finally, each 
Fourier amplitude will be expanded in a power series in e, with each 
coefficient of the power series a function of n. These functions of n 
are then found to be solutions of ordinary differential equations. 
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First, assume that 

♦ «,n,T> . I [2<„ ♦ 1 ) - tj . 2 „ tl 5 2n+l Wn.T). 
n»u 

Substituting (6) and (17) into (14) and equating the coefficient 
successive powers of 5 on the left and right sides of equation 
partial differential equations for the (F.(h,T)} are found. F 


51a.* gjSL.p<t,«,.> 

an 3 *3n 2 \3n / 3t3n 


Q(t;€,o), 


where 


P(t :c ,o ) - eo sin or , 


and 


Q(t ;e ,o) - 1 + i e 2 ♦ 2e cos ot + i e 2 cos 2 ot, 
2 2 

Next, define the operators 


3 3 F. 3 z F. 3F 3F . 3 2 Fj 

D(F.,F i; j) 5 —l + Fi— 2- (j + 1) —1—1 ♦ j F. 


3rr 


3n- 


3h 3n 3n 2 J 


and 


G(F. ,F . ) = — - 


3 F. 3 F , 3 2 F , 


F. . 


1 J 3n 3n 3n 2 J 


( 17 ) 

of 

(14), the 
satisfies 


(18) 


(19a) 


(19b) 


3 2 F , 

i (20) 

3x3n 


( 21 ) 
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,Fi}3) - i P - q, (22) 

4 

D(F5.Fl ,5) » i P - (1 4- i C 5 ) Q 4- 8C 5 G(F 3 ,F 3 ), (23) 

6 2 

D(F^ ,Fj ;7) -ip- (1 -f C7) Q ♦ 9C 7 G(Fs,F 3 ) 4- 15 C 7 G(F 3 ,Fs), (24) 

8 

D(F 9 ,F l ;9) - — P - (L ♦ C 4- 1 C ) Q + ic [48 G(F 7 ,F 3 ) 

LO 9,1 2 9,2 5 9,1 

4- 112G(F 3 ,F 7 )] 4- 18 C G(F S ,F S ), (25) 

9,2 3 3 

D(F n ,Fi ; 11) - i_ P - (1 + C + C ) Q 4- C [30G(F 3 ,F 9 ) 

12 11, 1 11,2 11,1 

4- L0G(F ,F )] 4- C [28G(F 5 ,F 7 ) 4- 20G(F 7 ,F s )J, (26) 

9 3 11,2 

with 


C S ■ /( fl l *5 ) 1 

(27a) 

" «3 «5 /(ai «7 )» 

(27b) 

°9,1 " «3»7/(aia9 ), Cg t 2 ■ a^/^ag), 

(27 c ,d) 

C 1 1,1 * *3 *9 / ( a l an), C llf2 ■ asa^a^ \) . 

(27e,f) 
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The boundary conditions for the (^n-fi^ * re 


F 2n.i<°' T) . <9F 2on /an) n . 0 - 0 


( 28 ) 


and 


[2(„ ♦ D-«J(»W»>* 1 ♦ e cos ox, 


(29) 


as . n + • . 


We have previously found the solution for Fi up to, and including, 

terras in e 3 (Grosch and Salwen, 1982; see Attachment) . Although the 

boundary layer equations are being used in this study, it should be noted 

that 5 Fj (n ,t ) is a solution of the Navier-Stokes equations for oscillating 

stagnation point flow. It is easily seen that the form of the solution for 

each of the F„ , is the same as that for Fi . Therefore, we have for 
2n+l 1 

j - 1,3 11 


Fj(n ,t) 


I 

n*0 


+ J, <**> 


+ 8 j,2m,k (,l) 8in (WT,! • 


(30) 


It is quite straightforward, but somewhat tedious, to obtain the ordi- 
nary differential equations and boundary conditions for the f's and g’ s. 
These differential equations are given in Appendix A. 

One of the major problems associated with series solutions of this type 
is that of determining the radius of convergence. On the basis of the 
results from the study of oscillating stagnation point flow (Grosch and 
Salwen, 1982), it can be concluded that the power series in e is rapidly 
convergent for e < 1, except perhaps at very low frequencies. If the fre- 
quency of oscillation is low, then a quasisteady approximation is quite 
accurate. Because of the rapid convergence of the series in e, for e 
small, only the terms through e 2 have been calculated. 
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A much note difficult problem is estimating the radius of convergence 
of the power series in It appears that there are no general results for 
determining the radius of convergence of the steady Blasius series, Van 
Dyke (1957) has shown that the steady Blasius series for a parabolic 
cylinder diverges for x/r Q > 0.62, where r Q is the radius of curvature of 
the parabola et the nose, The results given by Schlichting (1979) for the 
steady Blasius series for a sphere of radius r Q suggest that this series 

is convergent up to at least x/r Q " tt/ 2, It seems reasonable to assume that 
the time dependent series solution given here is convergent for € ® 7/1 & 1. 

RESULTS 


Numerical integration of the equations given in Appendix A yield the 
{ f, . ,g, . . ) , Substituting (30) and (17) into (13) and using (12) 

J | All | K ]j *nK 

gives, for the velocity components (u,v)» 


2n+l - , 

+ " S oJ a 2n+L 5 JJ C f 2n+l,2m,0 


'o ■ l {[ 2< " 

n-0 V 

+ Jj c2 "* ktf 2n + l > 2n,,k co “ (lM ’ T> + *2n + l,2m,k aio(ltOT)) ]} - < 31 > 


.-1/2 " | r . , . ^2n 


v/U^ - - (aj Re) £ {[ 2(n ♦ 1) - « o J(2n ♦ 1) a, nAl 5 


n«0 


2n+l 


I „ + I ‘"'U 


2mfk 


n*0 


?n+l,2m,0 v 2n+l,2ra,k 


cos (kot) + 


+ ®2n+l,2m,k ain 


(32) 


The velocity field depends on the body shape, the {aj} , the 
Reynolds number, Re, and the dimensionless frequency, o. it should be 
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noted that the effective Reynold* number and frequency depend, to tone 
extent, on the body shape because the velocity scale is *iU Q , see 

equations (8), (10), (11), and (13). It is clear from (31) that, as usual, 
the tangential component of the velocity, u, has no explicit dependence on 
the Reynolds number which only affects the scaling in the normal (n) 
direction. Again, as usual, the normal coiqponent, v, is proportional to 
d e ~l/2 £ n addition to the scaled dependence of the normal coordinate. 

Thrcugh the term* in e 2 , 


2 

u/U " ) + e u 

O 0,0 2,0 


(C ,n ) 




CO* OT + 




x sin ot] + c [u^ cos 2ot + u^ sin 2ot] , 

0,2 0,2 ' 


(33) 


v/U ■- (a.Re)” l/2 {v (5 ,n ) h v, (e ,n ) + e[v^ l J(5,h) cos or 

O 1 0|0 Z,0 0>i 


+ v< 2 } sin otj + e [v^(5,n) cos 2ox + v!? 2 J($,n) 

OpJL OpsC v 


(34) 


sin 3 jt] } 


with, 


u. 


J 0 I 2< " * l > - 5 oJ *2:i+l 


(35 a) 


u ( M - l [ 2 ( n 4- 1) - 5 ] a. . 5 Zn+l f* , .(n), 

n l “ V > ' nn J *) ni 1 Onil /» l' 1 '* 


2n+i 


o,J 


n"0 


on J 2n+l 2n+i,o,j 


(35b) 


( 2 ) 
u' . 
°>] 


H2(„ ♦ 1) - 6J . 2n+l E 2n+ ^„ +li0ij (n), 


(35c) 
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j,o “ L t2< “ * l) ' 'on' (2n * « *2n+l ^-l. J,o<"> • <”<*> 


J*0 


n«0 


,<*> 

o»J 


!* I [«n + l) -«J (2n + 1) ., ni , „ .(n). (35.) 


n*0 


on 


2n+l ’ 2n+l,o, j 


“j " t 2 <" + l > ' *.J < 2 " * l > *2n+l £ ^°®2n+l ,o, j £!1 ^ ’ < 35£ > 

n*o 


The superscript (1) indicate* the component of the velocity in the 
boundary layer which is in phase with the free stream oscillation, and the 
superscript (2) denotes the out of phase component. 

Figure l is a sketch of a typical blunt body and the coordinate system. 
The tangential and normal velocity components are u and v. The Cartesian 
coordinate system (x,y) is also shown in figure 1. The velocity compo- 
nents (u,v), parallel and perpendicular to the body axis are related to 
(u,v) by 


U ■ 

u cos 6 

- v sin 9 , 

(36 a) 

<> 

■ 

u sin 9 

+ v cos 6 . 

(36b) 


There is an inherent difficulty in presenting the results of the calcu- 
lations in such a way as to make them intelligible. The reason for this 
difficulty is that the character of the flow, in the boundary layer varia- 
bles, changes drastically with position. In front of the body the flow is 
primarily directed towards the body, so u is very small and v dominates, 
in fact for n large, v * - U 0 , As the flow moves around the body V 
decreases and u increases so that near the midpoint u dominates. For 
this reason it seems most reasonable to present the results in the (x,y) 
coordinate system. In particular, in terms of u, the component of the 
velocity is parallel to the axis of the body. Note, however, that as y is 
varied at fixed x both 5 and n vary. Presentation of the results in 
terms of u also facilitates comparison of the results of these calcula- 
tions with the experimental results of Kachanov, Kozlov, and Levchenko 
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(1978) although, aa will ba discussed below, the correspondence between the 
conditiona of the experiment and of thia calculation ia not exact, 

Detailed ealculationa have been carried out for the flow peat two types 
of bodies; elliptic and parabolic cylindera, 

% 3.1 Elliptic Cylindera 

On the aurface of an elliptic cylinder the tangential component of the 
potential flow ia 

H<5) - [l ♦ (b/a)] an(( ;k) , (37) 

with a the semi-major axis, b the semi-minor axis, 
an the Jacobian elliptic function, and 

k? - 1 - (b/a) 2 (38) 

The length scale ia the semi-major axis, a. 'Hie coefficients in the power 
aeries expansion for H(5), the (* 2 n+l^ * r * ** v * n * n 1» (Cayley, 

1895), Xhe radius of convergence of the Blasius aeries for an elliptic 
cylinder ia, apparently, not known. However, the ealculationa appear to 
converge up to at least the midpoint on the body, that is to 

S - E(k), (39) 

where E ia ihe complete elliptic integral of the second kind. The 
calculations were ended at this value of £ because we are primarily 
interested in the flow on the nose. 

Calculations of the velocity field have been carried out for a nunber 
of bodies with different slenderness ratios (a/b) and frequencies, o. 

Results are given for three bodies; a moderately thick ellipse with a/b ■ 

5; a thin ellipse, a/b ■ 10; and a very thin ellipse, a/b ■ 25, see 
Table 2, For each of these valuer/ of a/b calculations were carried out 
for a low frequency o *ir/4, two intermediate frequencies, o • ir and 
3n , and a high frequency, <r ■ 6* , In all cases the Reynolds number was 
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taken to b a 100.0, The raaulta at other, particularly higher, Reynolds 
numbers appear to hava tha aama form whan distances ara raacalad to tha 
appropriate boundary layer thickness. 

Contours of tha zero order steady flow, u q ara shown in Figures 2, 

3, and 4. Tha format of these figures, and all of tha others is tha sane. , 
The body shape is plotted in tha lower right hand corner of tha figure and 
the flow is from left to right. For each body the semi-major axis is the 
same length. The body shape has not been distorted, however tha region 
outside of the body has been stretched in order to show the details of the 
flow in the boundary layer. The coordinate system outside the body is the 
(5 ,ri ) system which is shown in Figure l. The distance 5*1 corresponds 
to four of the major units of the scale shown on the edges of Figure 2 and 
all others. The distance n - 1, on the other hand, corresponds to one of 
these major units. Thus, for example, the coordinates of the point in the 
lower left hand corner of Figure 2 are 5*0 and n ■ 4, and those of the 
upper right hand corner of Figure 2 are 5 * E(/0. 96) * 1.05050, and n ■ 
4,2. Note, however, that although the (5,n) coordinate system is used, 
these are contour plots of u, the velocity component parallel to axis of 
the ellipse, 

The results shown in Figures 2, 3, and 4 show that the zero order in 

e, mean flow has the same general features for bodies of different 

slenderness. Ahead of the body there is a boundary layer merging into the 

decelarating potential flow. Near the midpoint of the body the flow is 

nearly parallel to the body and the boundary layer flow merges into the 

tangential potential flow. In between there is a kind of "bubble" where 

u is small. This extends somewhat in front of and above the forward 
°»° 

portion of the body. This region of small u q q occurs because both u q q 

and v are small and are at a substantial angle to the x axis. As the 

0,0 

fluid flows around the nose of the body u increase* but v decreases 
until near the middle of the body where it is 0(Re"^^’). Therefore, on the 
forward portion of the body, but away from the axis, u is increasing with 
5 but is still not large and v is decreasing. Finally, as the body 
becomes more slender, the steady boundary layer tends to become thinner. 
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Figures 5 through 16 are contour plotc of th« second order ataady flow, 
u 2 Q . Figuraa 5 to 8 ahow tha raaulta for an elliptic cylinder with 

a/b - 5 for a sequence of frequencies, a » w/4, if, 3w, and 6f , 

Figuraa 9 to 12 ahow similar raaulta for a/b • 10, while tha raaulta for 
a/b * 25 are given in Figuraa 13 to 16, Examination of these figuraa 
ravaala aome general patterns in tha flow. 

In all cases u 0 ia poaitive in tha lower portion of tha boundary 

Z | p 

layer and negative in the upper portion. At low frequenciea the maximum and 
minimum valuea of u,, q occur at the midpoint of the ellipse. Aa tha 

frequency increases the maximum tendi to ahift towarda the front of the body 
and the thickneaa of the regiona of both poaitive and negative aecondary 
mean flow tend to decrease. However, aa the ellipae becomea more a lender 
i.e., a/b increaaea, the region of appreciable u. tenda to be confined 

ZpO 

to the near midsection of the body. The abaolute magnitude of thia aecon- 
dary flow alao tenda to decreaae with increaaing o, Finally, at low 
frequenciea the max (u ) ia greater than |min(u )|, but at high fre- 

quenciea thia ia reveraed, These results show that the Reynolds stress 

generated by the nonlinear interaction of the fundamental oscillatory flow 

with itself cause the formation of a weak, secondary, shear within the main 

steady boundary layer. At low frequency the maximum shear occurs near the 

middle of the boundary layer. As the frequency increases, the position of 

maximum secondary shear shifts towards ihe boundary. 

Contour plots of the amplitude and phase of u ., the fundamental 

o, l 

oscillation at frequency 0, are given in Figures 17 through 28. Aa above, 
these show results for a/b ■ 5, 10, and 25 anti a ■ rr/4, ir, 3n, and 6w. 
The phase ia, of course, arbitrary to within an additive constant. It has 
been assumed that for both the fundamental and the first harmonic, the phase 
is zero at $ ■ 0, n ■ 4, 

In all cases the contour maps of the amplitude of u , Figures 17a 

1 > 

to 28a, bear an obvious generic relationship to the corresponding contour 
plots of u for the same values of a/b. Of course, the structure of 

the flow does change with the frequency of oscillation. The amplitude of 
the oscillation decreases in front of the body and a region of reduced flow 
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is present ebove and, generally, in front of the nose of the body. As the 
oscillations move up onto the body, the amplitude grows. In all cases a 
peak in the amplitude distribution appears in the middle to lower part of 
the boundary layer at the midpoint of the ellipse. For a fixed value of 
a/b the location of this maximum of u, moves towards the boundary as the 

frequency increases, reflecting the development of a Stokes shear layer at 

high frequencies. If the frequency, a, is held fixed, the location of the 

maximum of u, also moves deeper into the boundary layer 4s a/b 

l,o 

increases, i.e. as the body becomes more slender. This is shown in Table 3, 

in which the location of the maximum of u , is tabulated as a function 

o, I 

a and a/b. 

The contour plots of the phase of u Figures 17b to 28b, are all 

o, l 

quite similar. The region of deceleration of the potential flow and the 

development of the boundary layer are obvious, as is the thinning of the 

boundary layer as it develops into a Stokes layer at high frequencies. Just 

back of the nose, the lines of constant phase tend to be parallel to the 

body with a noticeable thickening at low frequency near the midpoint of the 

ellipse. As a increases the region over which the phase shift occurs 

tends to thin. In all cases this region lies deeper in the boundary layer 

than the position of maximum u .. Independently of the value of a/b, 

o, i 

the total phase shift appears, from the results given in Table 4, to be 
tending to ir/4 as a + °* . 


Finally, Figures 29 through 40 are contour plots of the amplitude and 
phase of u q ^ , the first harmonic of the free stream oscillation. Again 

these results show very substantial internal consistency, 


The boundary layer in front of the body is readily apparent. The 
region of nonzero |u | decreases as the frequency increases at constant 
a/b and as the slenderness, a/b, increases at constant frequency. The 
amplitude of u q ^ has a maximum in the lower part of the boundary layer. 
For fixed a/b the position of this maximum moves lower in the boundary layer 
and forward on the body as o increases. At constant frequency the height 
of the maximum remains nearly fixed while moving back toward the midpoint of 
the ellipse as a/b increases. 
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The phase shift of u 
containing the maximum of 


o,2 


is generally nearly constant over the region 
and then goes negative in the outer 


portion of the boundary layer. In this outer region the phase then 
increases towards zero. At low frequencies this outer phase shift occurs 
over a thin layer, but at higher frequencies it is fairly broad, 


The first harmonic of the free stream oscillation is always confined to 
the steady boundary layer and, as a increases, is confined to a thin 
Stokes layer adjacent to the body, As the fluid moves up onto the body, 
away from the nose, the amplitude increases but never becomes large. In 
fact, |u .] is never greater than about 5 * 10" 2 . The velocity % 

°* z . i! 

fluctuations due to harmonics of the free stream oscillation thus never 
exceed 5 x 10“ 2 e 2 . 

3.2 Parabolic Cylinders 


We next consider the flow past a parabolic cylinder 


y 2 * 4ax, (40) 

where a is the distance from the nose to the vortex. Choosing the length 
scale to be a, the tangential component of the potential flow is a power 
series in odd powers of 5 » x/a. The coefficients of this series are given 
in Table 5. The solution does not contain any geometric parameters, that 
is, apart from the scaling, is independent of a and is, therefore, s 
similarity solution for all parabolas. 

The radius of curvature at the nose of a parabola is 

r o * 2a. (41) 

Using Van Dyke's (1957) result that the Blasius series diverges for x/r 0 
> 0.62, the solution given here is expected to converge only for 5 < 1,24. 

The coefficient ai * 1/4 for a parabolic cylinder, so that 
o " 4(toj l/ v 0 ). Thus the "effective" frequency in the differential equa- 
tions for the f's and g's is four times the true frequency of oscilla- 
tion. Of course, a similar scaling applies to elliptic cylinders, but for 
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an ellipse a^ ■ (1 + b/a)“i , which is nearly unity for slender ellipses. 
Because of the factor of four in the value of o for a parabola, 
calculations were carried out for o - if/16, if/4, 3n/4, and 3 tt/ 2, so 

that results would be readily comparable to those obtained for the ellipses. 

Contours of u Q ^ , the 0(1) mean flow past a parabolic cylinder are 
given in Figure 41. The format of this, and all of the other figures 
showing the results of the calculations for a parabolic cylinder, is similar 
to those figures showing the results for elliptic cylinders. Hie body is 
shown in the lower right hand corner. Distances normal to the boundary are 
scaled so that one major unit of the scale on the border of the figure 
corresponds to one unit of , the boundary layer coordinate* Thus the 
point in the lower left hand corner of the figures has the coordinates 
5 * 0.0, n ■ 6.0, and the top of the parabola at the right hand side of the 
figures has the coordinates 5 * 0.90, n » 0.0. 

A parabola is, in some ways, a rather peculiar blunt body in that it is 
not a closed body, so that there is no trailing edge, and that the thickness 
of the body is a monotonically increasing function of the distance from the 
stagnation point. Therefore, the farther back from the nose, the thicker is 
the region of potential flow which is substantially influenced by the body. 
This effect is apparent in Figure 41. There is a region of reduced flow in 
front of the parabola due to the deceleration of the potential flow and the 
formation of the boundary layer. This region extends far above the axis of 
the body because of the blocking effect of the parabola on the potential 
flow ahead of the body. In the region in front of the nose the contours of 
u are nearly equally spaced. Just above the parabola, these contours 
tend to line up with the body and the distance between them becomes 
increasingly smaller as the flow moves onto the parabola and the quasi- 
parallel boundary layer flow develops. 

. Contours of the steady, second order, streaming flow, u 9 , for 
Acr "it /4, it, 3it , and 6 tt , are shown in Figures 42 through 45. These 
plots show features which are similar to those seen in the contour plots of 
N « for the flow P ast an ellipse, Figures 5 through 16. In both sets, 

4,0 
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tig 0 is positive in the lower portion of the boundary layer and negative in 

che upper part. The thickness of the secondary boundary layer decreases 
aa the frequency increases and there is a slight tendency for the position 
of the relative maximum and minimum of u. to shift forward towards the 

4|0 

stagnation point as 0 increases. However, as o increases the maxlu. I 

2,0 

decreases significantly. 


Figures 46 to 49 are contour plots of the amplitude and phase of u q 

the fundamental oscillation in the free stream flow past the parabola. The 

contours of the amplitude of u , show a general similarity to those of 

o, l 

steady flow, u 


0,0 


However, as o increases a Stokes layer develops on 


the parabolic cylinder and moves forward towards the stagnation point. Away 
from the stagnation point, the contours of u , are similar to those of 

. °i 1 

for an ellipse if the surfaces are rotated so as to be locally 


o,l* 

parallel. The local maximum of u q seen in the flow past an ellipse is 

not, however, seen for the flow past the parabola. The contour plots of the 

phase of u ., also show the development of a Stokes shear layer. The 

o, l 

total magnitudes of the phase shift for u , in Figures 46b through 49b 

o, l 

are quite close to those for the corresponding cases for the elliptic 
cyi inders. 


Contour plots of the amplitude and phase of 



Che first harmonic 


of the oscillation in the free stream speed, are shown in Figures 50 to 53. 

The results shown in these figures illustrate the development of a second 

order Stokes shear layer within the steady boundary layer as the frequency 

increases. The raaxju | increases with £, all o, at least up to 

o ,2 


the point £ “0.9. It cannot be determined whether or not there is a 
localized maximum in the amplitude, as was found for the elliptic cylinders, 
or whether the amplitude reaches a plateau. The problem is that, in order 
to decide this question, the calculations would have to be extended sub- 
stantially beyond £ *0.9, but, as mentioned above, the series expansion 
probably is divergent for K > 1.24, and probably converges very slowly for 
£ close to 1.24. 
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DISCUSSION AND CONCLUSIONS 


Perhaps the most general conclusion which can be drawn from the results 
of this study is that in the region of the nose of a symmetric, two dimen- 
sional blunt body at zero angle of attack, the steady plus oscillating flow 
is very similar for a wide class of body shapes, this conclusioon has been 
shown to be true for elliptic cylinders with a/b < 25, and for the parabo- 
lic cylinder. Additional calculations, not reported here, were carried out 
for elliptic cylinders with values of a/b up to 100, with results which 
are very similar to those reported here. In all cases, the flow field in 
the nose region of a two dimensional blunt body is generic to that of the 
flow in the neighborhood of the stagnation point on a plane wall. 

A general picture of the flow on the forward portion of a blunt body, 
due to a steady plus osciLlating free stream, can be sketched. Forward of 
the body there is a region of decelerating potential flow which merges into 
the viscous stagnation point boundary layer. On the stagnation stream Line 
the boundary layer thickness is 2.3&t (Re)” 1 ^ 2 , with t the scale length 
of the body. If one considers u, the velocity component parallel to the 
body axis, there is a region of small u, above and ahead of the nose, 
where the velocity is reduced and is at a significant angle to the body 
axis. On the body itself, the boundary layer changes slowly with distance 
from the leading edge. The description, so far, applies to u , the 0(1) 

O |0 

steady flow. 

If one next turns to u , the fundamental oscillating component of 

o,l 

the flow, the above description is, with some additions and emendations, 
valid. These are due to the changes with frequency in the oscillating flow. 

At low frequencies, mi /a L U < — , the flow is essentially quasi-steady; 

° 4 

that is, it is the steady flow scaled by the instantaneous free stream 
speed. As the frequency increases, the oscillating boundary layer develops 
into a Stokes shear layer. This change occurs along, at least, the entire 
forward portion of the body, independent of the position on the body. This 

is not completely obvious when u , is examined, but it is obvious if 

O' , /* 
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studied. As shown by Grosch & Salwen (1982), Che development of a Stokes 
layer at high frequency occurs at, and in the neighborhood of, a stagnation 
point. 

This is in marked contrast to the results of Ackerberg and Phillips 
(1972), who concluded that the Stokes layer only develops far downstream of 
the leading edge of a flat plate. This may be due to the fact that' there is 
a singularity at the leading edge of the (nonphysical) infinitely thin flat 
plate, or it may be due to the scaling assumptions of Ackerberg and 
Phillips, which, in effect, equate low frequencies with small distance from 
the leading edge and large distance from the leading edge with high frequen- 
cies. In the calculations reported here, the distance from the stagnation 
point and the frequency are independent. 

The oscillations in the boundary layer, which exist at all frequencies 
and at all positions on a blunt body, could be interpreted as an oscillation 
in the steady boundary layer thickness. In the free stream, the mean (time 
averaged) speed is U q , and the boundary layer thickness is conventionally 
defined as the height in the boundary layer, 5, where u(5) ■ 0.99 U . As 
u increases and decreases at fixed n, because of the oscillation, the 
instantaneous value of 6 would appear to increase and decrease. It is 
quite easy to estimate the amplitude of the apparent oscillation in 6, if 
the velocity oscillations in the free stream are small. Let 

5 - 5 q + 6i(t), (42) 

with S Q the time averaged boundary layer thickness, and 5i the oscillat- 
ing component of the boundary layer thickness. By definition, 

<u(5)> - 0.09 <U> - 0.99 U , (43) 

o 


where <> is a time average. Then, if e « 1, it is easy to see that, to 
lowest order in e, 


3u 

3 y y * 6 


e /( 

' ' rw 'Vs# 


cos (wt + 0 ) , 


0 J 


(44) 
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where, as above, u ia the 0(i) steady boundary layer flow component 


0,0 


parallel to the boundary, and 9 is the phase of m q evaluated at 

~ t * 

y ■ s 0 - 

In the neighborhood of the stagnation point on any blunt, two dimen- 
sional body, 6, scaled by l, is 


«o - 2.38 (Re)“ 1/2 * 


(45) 


and 


5u 

6 y y ■ 5 


2.6 x 10- 2 . 


(46) 


So that, in the vicinity of the stagnation point, |$i|i scaled by l, is 


|6 ! | *• 38.5 e( ai Re)" 1/2 


(47) 


and 


j /6 o - 16.2 e( ai )- 1/2 . (^8) 

At low frequencies the flow is very nearly quasi-steady, but at high 

frequencies, there can be some distortion of the velocity profile. If o 

is large, u is essentially constant outside of the Stokes layer which 
o,l 

has a thickness of 0(6 o -1/2 ). Inside the layer, u . can cause some 

o o,i 

distortion of u , while outside the Stokes layer u . only adds a 

o,o o,I 

constant, in n , time varying increment to u Qj0 . 

It may be useful, at this point, to give an idea of what is a low 
frequency and what is a high frequency. From equation (11), 

a m ml / (*i U 0 ) . 

If, by way of illustration, l ■ Ira. , ai ■ 1, and Uo ■ 10 m. /sec. , then 
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1> sl» 


0 - *f/5 (49) 

where £ is the frequency in Hertz. Thus for 

o .$ ff /4 , 

the low frequency region, 


f * 5/4 Hertz. 


While, for high frequencies, 
o Z 6jt , f Z 30 Hertz. 

In summary, for this example, frequencies below 1 Hz. are low frequencies, 
those greater than 30 Hz. are high frequencies, and the range of 1 to 30 Hz. 
is the intermediate range of frequencies. 

Considering next the terms of 0(e 2 ), thare is both a second order 
steady streaming flow, u , and u , the first harmonic of the free 

4)0 O) 4 

stream oscillation. The steady streaming flow extends throughout the 0( L) 
boundary layer at low frequencies, but at high frequencies, it is largely, 
but not completely, confined to a Stokes layer. However, at high 

frequencies, there is a small portion of which decays to zero in an 

* 

outer layer, with a thickness of the same order as that of the 0(1) mean 
boundary layer. 

The magnitude of this flow, |u |, increases with £, and for an 

elliptic cylinder has a maximum just in front of the midpoint of the 
ellipse. Such a maximum may, or may not, occur on a parabolic cylinder. 
Because of the limited range of convergence of the series expansion, the 
calculation cannot be reliably extended far enough to decide. 

^ The first harmonic, u q has a structure which is partially similar 
to that of the secondary streaming flow, u 2 D » and partially, to the 
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fundamental, u q At low frequencies it extends throughout the steady 

boundary layer and at high frequencies is confined to a Stokes layer. It 
also develops a maximum amplitude some distance from the stagnation point. 

Both the first harmonic and the secondary streaming flow are present at 
all 5 , even at the stagnation point. Away from the stagnation point they 
are somewhat modified, but have the same character as in the neighborhood of 
the stagnation point. In particular, the estimate (Grosch & Salwen, 1982) 
that they are bounded by (e/0) 2 seems to be valid over an appreciable 
range of 5, in fact, up to the midpoint on elliptic cylinders and over the 
range of a convergence of the series for flow past a parabolic cylinder. 

The only experiment with which the results of these calculations can, 
apparently, be compared is that of Kachanov, Kozlov, and Levchenko (1978), 
They studied the flow past a flat plate with a nose consisting of two con- 
jugate ellipses with a semi-major axis, a ■ 50 mm. and semi-minor axis, b ■ 

2 mm. on the working side and a - 128 mm. and b ■ 8 mm. on the other side. 
Figure 54 shows the cross section of this body from the nose to a distance 
of 20 mm. back from the nose. 

In Figure 2 of their paper, Kachanov et, al, give cpntours of (in the 

notation used here) u in the region from about 10 mm. in front of the 

O } o 

nose to about 15 mm. behind it, and from the axis of the body to 6 mm. 
above. Figure 55 is an enlarged copy, supplied by Dr. Levchenko, of Figure 
2 of their paper. In this figure the contours are labeled with the speed in 
units of meters/second. 

It is obvious that there is a general, qualitative agreement between 
the theoretical results shown in Figures 2, 3, and 4 for the mean flow and 
the experimental results of Kachanov, et, al. shown in Figure 55. The major 
difference between the theoretical and experimental results is the closed 
contour labeled 6.0 in Figure 55; nothing like this is seen in the results 
of the calculations. The experimental result is rather curious. Either 
there is an absolute maximum within this contour, a "peak", or a local mini- 
mum lies within the contour, so that the contour is the "lip of a volcano." 
In either case the flow speed, as shown in Figure 55, is a maximum outside 
of the boundary layer. 
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It is well known that, for a potential flow, the maximum velocity 
occurs on the boundary, The potential flow past an ellipse is a classic 
problem (Milne-Thomson, 1955). Calculations of the potential flow for an 
ellipse with a/b ■ 25 confirm that the maximum of u occurs on the 

Of 0 

boundary at the midpoint of the ellipse. Taking into account the existence 
of the viscous boundary layer, the maximum speed should occur at the top of 
the boundary layer, How is it then possible to account for the existence of 
a maximum in the speed in what is, apparently, a region of potential flow? 

There seems to be three possibilities, The first, and most obvious, is 
experimental error. This does not appear to be likely because it would 
require systematic errors of the order of 102 to 152. A second possibility 
is that the stagnation point is not on the axis of the body (see Figure 54). 
If, as seems likely, the stagnation point lies below the axis, then the flow 
must first move through a region of adverse pressure gradient until it gets 
around the nose. Then, for a short distance, the pressure gradient is 
favorable, and then falls rapidly to almost zero just behind the nose. It 
is possible that there is a local flow separation and reattachment, a 
separation bubble, just behind the nose of the upper ellipse. This would 
lift the boundary layer in this region and might account for the closed 
contour of u in the results of Kachanov, Kozlov, and Levchenko. The 

O | o 

third possibility is that the free stream vorticity is not zero. If the 
free stream vorticity is non-zero, perhaps there is % shear in the free 
stream, then, depending on the distribution of free stremn vorticity, a 
maximum of the speed could occur virtually anywhere outside of the boundary 
layer. 

There is, of course, a final possibility. These theoretical argisnents 
and calculations may be wrong. It may be possible, for some reason, that a 
noir* rotational potential flow has a maximum away from the boundaries. Or, 
it may be that there is a viscous effect, not included in these calcula- 
tions, which causes a maximum in the speed outside of the viscous boundary 
layer. 

Kachanov, et. al. also give experimental results for an oscillatory 
past this body. The oscillatory flow is caused by a vibrating ribbon above 
the axis of the plate and in front of it. As Kachanov, et, al, indicate, 
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the periodic vortices shed by this ribbon passed above the plate and only 
the tails of the vortices impinged on the nose of the plate. The results of 
the calculations given here are, at best, in rough qualitative agreement 
with the experimental results of Kachanov, Kozlov, and Levchenko, In pert, 
their results show an absolute maximum in the amplitude of the oscillatory 
velocity in a small region just above the axis of the plate and just beyond 
the top of the steady boundary layer. They also show a rapid change in the 
oscillatory flow inmediately below the region of maximum which wea discussed 
above. As was discussed above, either the calculations reported here are in 
error or the experimental results reflect unnoticed effects. In fact, from 
the results given in Figure 1 of Kachanov, et. al, it appears that the 
oscillation has a nonzero shear in the free stream. In any event, this 
issue requires clarification and would seem to require further work, both on 
the experimental and theoretical aspects of this problem. 

Finally, if the theoretical approach used here and the results are 
correct, it would seem to have value beyond this particular study. The 
general character of the flow, in the nose region of the blunt bodies con- 
sidered here, is very similar to that of the flow in the immediate neighbor- 
hood of the stagnation point. As the flow moves around the nose of the 
body, no essentially new features appear. This suggests that the structure 
of the flow on the forward part of a blunt body due to more general classes 
of disturbances, such as vorticity waves, can be found by studying the flow 
due to these disturbances in the neighborhood of a stagnation point. It 
seems clear, at least, for two dimensional flows, that a linearized, in e, 
theory is sufficiently accurate. The results given here show that, for any 
ellipse and for any parabola, the second order terms are bounded by (e/o) 2 
for oil. It also should be noted that the steady flow in the neighbor- 
hood of a general three-dimensional stagnation point is known. This should 
facilitate the theoretical study of disturbances impinging on a body in the 
region of a stagnation point. These approximations, flow near a stagnation 
point, and linearization in c, should be a considerable simplification for 
both analytical and numerical studies of boundary layer receptivity. 
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APPENDIX A 


Differencial Equations and Boundary Conditions for the F'a and G's 

The function F^ (n ,t ) ia the solution of the oscillating stagnation 
point flow problem (Groach & Salwen, 1982; see attachment). The differen* 
tial equations and boundary conditions for the f’s and g'a are given 
there through terms e 3 and will not be repeated here. 

The differential equations and boundary conditions for the f's and 
g'a given here are those in the expansion of the {F^(n,T)} for 
J ■ 3,5, ... ,11 and e 11 , with n £ 2. First, we define an operator L by 





<j> W 



(Al) 


where t , (ti) is the 0( L) function in Fi , and is the Hiemenz function, 

and primes denote differentiation with respect to n. Note that in (Al) the 
coefficients are j and j + 1 and that j is the first index of the f 
or g upon which L operates. Next we define a second operator S by 

S(f(n), g(n);a,M) 5 afg" - $f'g' + Yf" g, <A2) 


where f and g are any functions of n • primes, again, denote differen- 
tiation and a, 0, and Y are arbitrary constants. 

Then for the terms proportional to 5 3 we have, at 0(1) 

Lf, ■ - l (A3a) 

3,o,o 
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f, r. r* « (0) " °» 

3,0,0 3,o,o 


i 


3,0,0 S' 


+ 7 - aa n + 


At 0(e) we have 


W 3,o,l ‘ 'S.o.l ' * l - S( ‘l,o,l ( , 3,o ) o il -' ,,3)> 


l 


Lg 3,o,l + 3 f 3 , o , 1 “7° ' S(g l,o,l* f 3,o,o J l|4,3) » 


f 3,o,l <0) " f 3,o,i (0) " g 3,0,l (0) " g 3,0,l <0) " °* 


*3,0,1 + 7- *3,0,1 + 0 “ " 


And at 0(e 2 ) 


Uf 3,2,o " " 7 " S(f l,2,o’ f 3,o,o’ l > 4 ’ 3) 


" I t S<f l,o,l ,f 3,o,l ;1,4>3) + S(g l,o,l’ 8 3,o,l 5l,4 » 3 ^ * 


f 3,2,o<°> " f i,2,o (0) “ 0 


f’ ♦ 0 as n ♦ 
3,2,o 


and 


Lf - 2 0 a' 
3,o,2 8 3 ,o,2 


‘ S(f l,o,2* f 3,o,o ;l ’ 4 ‘ 3) 

- “ [s(f l,o,l ,£ 3,o,l il, *’ 3) ' S( *l,o,l>*3,o,r l '*- 3) ]- 


(A3b) 

(A3c) 

(Ma) 

(Mb) 

(Me) 

(Md) 

(A5a) 

(A5b) 
( A5c) 

(A6a) 
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Lg 3,o,2 * 20 ^3 ,o,2 8<g l,o,2 ,f 3,o,o' ,l,4) ^ 

“ 7 t S ( f l,o,l ,g 3,o,l ;l,4> 3) + 8(g l,o,l» f 3,o,l {l » 4,3) l» (A6b) 


f 3,o,2< 0> " f 3,o,2 (0) 


g 3,o,2< 0) 


3,0,2 


( 0 ) - 0 , 


f 3,o,2* g 3,o,2 + 0 *• n + 


( A6c) 
<A6d) 


Next the f'a and g' s proportional to 5® satisfy, first to 
“5,0,0 < l *}°5> ‘ 8 ° s S(f 3.o.o- 
£ 5,o,o«» ‘ £ 5,o,o«» * 0 
£ 5,o,o"£ « 

Then at 0(e) 

Lf 5,o,l " ®*5,o,l * " 2<l * -j c 5 ) " *«l,o,l> t 5,o,oi l .‘.» 

- *°5 

“5,o,l ' o£ 5,o,l ■}* - s <*l,o ) l- f 5,o,oi l > 6 - 5 > 

-4C 5 S(f3 >0j0 .*3 >0>1 !l.2,U, 

£ 5,o,l< 0 > ’ £ 5,o,l«» ' «5,o,l< 0 > - *5,o,l< 0 > * 


£ 5,0,1 * *5,o,l + 0 " * * 


The equations and boundary conditions at 0(e 2 ) are: 

k £ 5,2,o'" " 4 <1 * - »< £ l,2,0. £ 5,o,oi‘*«.S) - \ 


2 2 * M*|V- 2 

* [ S ^ f l,o,l» f 5,o,l ;1 * 6 * 5 ^ + S ^ 8 l,o,l» 8 5,o,l ; 1,6 » 5 ^ 

- 8 C 5 {s(f 3>2|0 £ 3>0t0 ;l,2,l) 


0 ( 1 ) 

(A7a) 

(A7b) 

(A7c) 


(A8a) 


( A8b) 
(A8c) 
(A8d) 


i \ 
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+ \ t S(f 3,o,l' f 3,o,l 5 * S( *3,o,Vh t o>V 


v l *i]) 


f 5,2,o (0) " f 5,2,o (0) 


0; f 5,2,o * 0 w n * **’ 


(A9a) 

(A9b,c) 


Lf 5,o,2 ” n 5,o,2 " ” 7 (l * 7 C 5> S<f l,o,2 ,f 5,o,o ;1 » 6 » 5) 

”7 f S<£ l,o,l ,f 5,o,i ;l ‘ 6,5) " S< 8 i,o,i , 8 5,o,l Jl ’ 6,5) l 


- SC. 


i( S < f 3,o,o’ f 3,o,2‘ l ’ 2 ’ l) + 7 [ S(f 3,o,l» f 3,o,l ; b l ’¥ 


S( 8 3,o,r 8 3,o,r 2 ,l, 2^) 


( AlOa) 


Lg 5,o,2 4 2 of 5,o,2 " “ S<8 i,o,2 ,f 5,o,o ;1 » 6 ’ 5) "7 [ S(f l,o,l’ g 5,0,l J 
x 1,6,5) ♦ S(8 1 # 0 > 1 f 5 > 0 > 1 ;1,6,5)] - 8 C 5 {s(f 3 >o>o ,g 3 >o>2 ; 

x 1 , 2 , 1 ) + ‘y[ 8 ^ 3 | 0 > 1 .* 3 | 0 fl 4 »l»’j) + { s( * 3 ,o,l’. f 3 ,o,l J 
x 7 , 1 , 7 )]], (AlOb) 


fc AO) - f' 9 (0) « g, ,(0) « g' ,(0) - 0; 

5 , 0,2 5 , 0,2 5 , 0,2 5 , 0,2 


’ f 5,o,2 * *5,0,2 * 0 “ n * ”• 


(A10c,d) 
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The f'a and g f a which are proportional to C 7 aatiafy the 
equationa and boundary conditiona, at 0(1) 

Lf - l (1 ♦ C ) - C S(f ,f .9,24,15), 

/ ,o v o 7 7 J ,o , o 5,o,o 


f, (0) - fi (0) - 0, fl ♦ i aa n + 

7 A A 7 A A ' ' • 7 A A ’ 


• ■ \v / a. 

7,o,o 7,o,o 


7,o,o 


8 


At 0(c) the equationa and boundary conditiona are: 


Lf - Oa' 

7,o,l 8 7 ,o,1 


- 2(1 + C ) - S ( f f ; 1,8,7) 

7 1 , o , 1 7,0,1 


-C ? [s(f _ _ > f e n i J 9, 24, 15) + S(f ,f ; 

7 3,0,0 3,0,1 3,0,1 3,0,0 

X 9,24,15)], 

L *7,o,l + ‘ £ Ko,l-;°- S( *l.o,l- £ 7,o,oi l ’ 8 ’ 7 >- 

* C 7 [S<£ 3,o,o- g 5,o,l> 9 * 24 - 15 > + S<8 3.o,l- £ 5 ,o,o> 
x 9,24,15)], 

f_ . (0) - f ' 1 (0) ■ g 7 _ ! (0) ■ . (0) - 0, 

7,0,1 7,o,l 7,o,l 7,0,1 

£ },o,l*J- *7,0,1 ^ 0 
While at 0(e 2 ) we find 


Lf 


7,2,o M_ “ (1 + C 5 ) “ S(f l,2,o’ f 7,o,o ;1 * 8 ’ 7) “ “t S(f l,o,l» 


f 7,o,l ;1,8,7) + S(g l,o.l’ 8 7,o,l ;1,8,7) ^ 


(Alla) 

(Allb,c) 

(Al2a) 

( A12b) 

(Al2c> 

(A12d,e) 


33 



ORIGINAL PAGE IS 

OF POOR QUALITY 

- C 7 {S<f 3,2 > o > f 5 > o,o' 9 * 24 ' l5) + S(f 3,o,o* f 5,2,o> 9 ' 24 ’ l5) 
+ 7t S(f 3,o,l' f 5,o,i ;9 ' 24 > l5) + S(g 3,o,l» g 5,o,i ; 


x 9,24,15)]), 


(A13a) 


f_ . (0) - fl (0) - 0, fi . + 0 as n (A13b,c) 

7,2,o 7,2,0 7,2,0 


L£ 7,o, 2 - tog 7,o, 2 ' - f * V - S<£ 1,o, 2 ' £ 7,o,o> 1 > 8 - 7) 

■ J [ S<£ l,o,l’ f 7,o,r * S<g l,o,l ,g 7,o,l ! 

* 1,8,7)] - C 7 tS(f 3 o o ,! 5 o 2 ;9,24,15) ♦ »<£,,., ,£,^.1 

> 

x 9, 2 4,15) * g <*3,o,l-*5,o.l> 


X 9,24,15)]}, 


(Al4a) 


Lg 7,o, 2 + * ' S<g l,o, 2 £ 7,o,o il,8,7) “ I [ S(£ l,o,l ,g 7,o,r 

x 1,8,7) * *« g ll ., 1 . £ ,^. l '*. 8 . 7 > - C 7 lS(£ 3,o,o- g 5, 0 , 2 ! 

x 9,24,15) ♦ g ^ g 3 ,o, 2 ,£ 5 ,o,o’ 9>7 ^ ,£87 * ?t S<f 3 ,o,l ,g 5 ,0,1 : 


x », 2 4,15> ♦ g < g 3 , 0(1 .*j >0i ,l». 2 4,lS)]>. 

£ (0) - £1 (0) - g (0) - gi (0) - 0, 

f>0,l 7,0,2 7,0,2 7,0,2 


(A14b) 


( A14c) 


f ' + gl „ ♦ 0 as n 

7 ,o , 2 7,0,2 


(A14d) 
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For the terns proportional to £ 9 the 0(1) equations and boundary 
conditions are: 


If, 


1 , 0)0 “ “ <l + C 9,i * I C 9,2 ) " Iff C 9,l S(f 3,o,o’ £ 7,o,o’ 


l . 1, 


96,320,224) - l8C 9 t 2 S ^ f 5 ,o,o ; f 5,o,o ’T* ’ 


f Q <0) m f' Q n JO) - 0, f ' o ♦ 1/10 as H + 
9,o,o 9,o,o 9,o,o 


(A15a) 


(A15b,c) 


The 0(e) equations and boundary conditions are: 

L£ 9,o, 1 " N.o.l ' _2(1 * C 9,l * Vi' ‘ S<f l,o,l ,f 9,o,o :l ' 10,,) 
' TO °9,J S(£ 3,o,o ,£ 7,o,l l96,320 ' 224> * S(£ 3,0,l ,£ 7,o,o i 


(Al6a) 


* 96,320,224)] - 1% j 3 ; 8 <« 5 ,o.o' * 5 , 0 , 1 ! l ’ 2 ’ l) 

L8 9,o,l * ° £ 9,o,1 " W 0 " S<8 l,o,l’ £ 9,o,o il>10 ’ 9) ’ Iff C 9 ,1 

* ^ 8££ 3,o,o ,8 7,o, l’ 94 ,320, 22 4) ♦ *<«,, o ,i* 7 , Ol .»»‘^0,22*)] 


- 18C. „S(f ,g, 1,2,1), 

9,2 5,o,o 5,o,l 


£ 9,o,l (0) - £ 9,o,l {0) - 8 9,o,l (0) - 8 ;, O , 1 (0) -°' 


f 9,o,l* TO* *9.0,1 * ° 08 n> "' 


The 0(e 2 ) equations and boundary conditions are: 


(Al6b) 
( A16c) 
(Al6d,e) 


Lf 


9,2,o 2 


- - i(i + c 9>1 + i c 9>2 ) - s(f 1 > 2 > 0 f 9>0>0 ; 
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» >.10.*)] -Tts(f 1>Oil .« wi l,l 0> 9)*,( gl|O(llg , iO>li 

« 1.10,9)] * S((, >o o . 

* £ 7 , 2 , O !96 . 320 " 224) *?ts(( 3 >o> 1 ,f 7jO l i 96 , 320 , 224 ) 

* S(g 3,o,l ,8 7,o,l !96 > 320 . 224) ]} ' l 8C 9,2{ S<f 5,2,o ,£ 5 > o I o i 

xi. 2.0 4[«f Si0iI f St0>I , * °<',, 0 ,rh.o,r 

X !■ 1 .yl ] ) (A17.) 


f 9.2,o (0) * f 9,2,o <0 > ' °" f 9,2.o * 0 “ " * *• 


(Al7b,c) 


Lf 


9,0,2 


- 2< S,o, 2 ' -T ( ‘ + °9.1 *T C 9,2> - S(f llO ,2. f 9.o.o = l.l 0 . 9 > 
' T [ S<f l,o,l- f 9,o,O l . l0 - 9) * S(8 l. O ,l- 8 9,o,O l . 10 . 9 >] 
-TC C 9,l{ S(f 3,o,o £ 7,o,2l 96 . 32 °. 224 > * S<f 3,o,2. 

* f 7 ,o,o’ 9 £, 32 0 ,224 l + yt S(£ 3 .o,l- £ 7 .o.O 96 . 320 . n4 > 

' S<8 3,o,l' 8 7,o,l i96 . 320 > 224> ]} ' 18C 9,2{ S<£ 5,o,o' f 5,o,2 ! 

x 1 , 2,0 4t S < £ 5 ,„,l £ 5. 0 .O7.1.y>-« 8 5 , 0 . 1 . 8 5,0,O 


(A18a) 


Lg 9,o,2 + 2 ° f 9,o,2 


" S(8 l,o,2» f 9,o,o ;l * 10 ‘ 9) ’ 1 ^ S(f l ,o,l * 8 9,o,l * 
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Si 

1,10,9) + 

s(g i,o,r f 9,o,i ;l,l °’ 9 ^ " 

1 

10 

C 9, 1 S(f 3,o,o’ 

X 

8 7,o,2 ;96 » 

320,224) * S(8 3f0>2 ,f 7>0|0 

; 96 , 

,320,224) + j 

X 

t S(f 3,o,l* 

g ? ^96,320,224) ♦ S(g 3> 

0,1* 

f 7,o,l> 96 ’ 320 ’ 

X 

224)]} - l 8 S >2 {s<f 5i0 _ 0 ,g 5i0i2 il,2, 

l) ■* 

’ S(f 5,o,l’ 8 5,o,t ! 

X 

1,2.1)}. 



(A18b) 

f 9,o,2«» ' 

' f 9,o.2<°> 

■ 8 9,0,2 (0) * *9,o,2 (0 > - 

0, 

(A18c) 

f* + 

9,o,2 g 9 

',0,2 

as n <*. 


( Aldd) 


Finally, for the terras proportional to C 11 the differential equations 
and boundary conditions are, at 0(1): 


Lf 


11,0,0 - - (1 * °11,1 + °ll,2 ) - C ll,l S(f 3,o,o- f 9,o,o: 10 - W - 30) 


- °U,2 S(f 5,o,o- f 7,o,o: 20 - 48 ' 28) 


(A19a) 


f ll „ „<°> ’ f ii „ „«» ■ 0, 

1 l ) O |0 11)0,0 


11 , 0,0 12 


as q + «. 


(A19b) 

(A19c) 


The 0(e) differential equations and boundary conditions are: 


Lf l 1 n 1 “ i i * " 2(1 + G,, , + C.. J - S(f. ..f., ; 

^ i ^ ^ 11,2 i,o,l 1 1 1 o ) o 


1,12,11) -C U I [ S (f 3 _ OiO ,£ 9iO|li 10, 40.30) 


-K s (f ,f ; 10 ,4°, 30 )j - c:.. _ [s(f_ ,f_ : 

3,0,1 9,o,l J 11,2 L 5,o, o’ 7,o,l 
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* * 


* 20,48,28)) ♦ S(f 5 ol ,f 7 O j20,48,28)], <A20a) 


L8 U,o,l + 5 f ll,o,l " 12°' " S(8 l,o,l ,f ll,o,o ;l ' 12 » 


il> 


' C ll,l^ S(f 3,o,o ,8 9,o,l ;l0,A0,30) + S(8 3,o,l ,f 9,o,o ; 
* 10,40,30)] - C U(2 [s(« 5 o>o ,g 7iO(li 20,48,28) 

+ S<8 5,o,l- f 7,o,o‘ 2 °'* 8 ' 28 >J- 


.( A20b) 




f n « i<°> " f !i „ i (0) " 8|1 , <°> " 8{, ,(0) - 0, 

ll,o, l lL,o,l ll,o,l ii,o,l 


* IT- g h,o,l 


0 as n , 


( A20c) 


(A20d,e) 


And, lastly, the 0(e 2 ) equations and boundary conditions are: 

ll,2,o " ‘ I (l + C ll,l * C ll,2 ) " S(f l,o,2 ,f ll,o,o 5l,12>ll) 

“ 2 ^ S<f l,o,l ,f ll,o,l ;l,l2,ll) + S(g l,o,l ,8 ll.,o,l i 
x 1,12,11)] - C Uil {s<£ 3>2(O .f, iO|o! 10, 40,30) ♦ S (£ 3i0 o , 

x f,_ 2|0 l 10,40,30) + -j t 8 ^ 3 , o,l *^ 9 , 0,1 *10,40,30) 

* S(8 3, O ,l- 8 9,o,l' 10 - W > 30) ]>- 0 U,2{ S(f 5,2,o- f 7, O ,o i 
x 20,48,28) * S(f 5>O(O ,£ 7j2ioi 20.48,28) * ‘ [s(f 5>ol , 

x £ 7>o>l i20,48.28) * S( g5Ojl ,g 7tOil ;20,48.28)]}, (A 21 .) 


*11,2, o ( °> ‘ f il,2,o (0) * °* 


(A2 lb) 


38 


ORIGINAL PAGE IS 
QF POOR QUALITY 


f * 

c lt,2,o 


o«§n* 




(A21c) 

' 2og ii,»,2 

.-ia 

+ C U,1 * 

C ll.2> 

" S(f l,o,2* f U,o,o ; 

x 1,12,11) 


,o,l' f U,c 

i,l !l ’ 12 

,10 - s(g lol , 

" *U,o,l il 

,12,11)] 

- C.. . S(f, , 

11,1 3,o,o’ 

f 9,o,2* l0 * 40 * 30) 

* S(f 3,o,2' 

f Q ;io 

9,o,o’ 

,40,30) + 

T t s(f 3 

,o,l* f 9,o,l* l0 * 40 * 30) 

* S(8 3,o,1> 

g 9,o,l ;10 

,40,30)] 

’ °ll,2 

S(f, . A) 

5 j o 1 0 7,0,2 

x 20,48,28) 

* S(f 5,o 

,2’ f 7,o,o : 

20,48,28) [ s < f 5|0|1 . 

X f 7,o,l* 20 

,48,28) - 

S<g 5,o,l- 

*7,o,l* 

20,48,28)] , (A22«) 

' 2af U,o,2 

■ ~ S(g. 

1 » 

o,2’ f ll,o, 

o’ 1 * 12 * 

11) 

"7 [s(f i,o 

,i’ g ll,o, 

*1*12,11) 

+ s(g t 

,0,1’ 

X f U,o,l ;l 

,12,10] 

■ °n,i! s< 

f 3,o,o'«9 > o > 2' l0 > 40 - 3< » 

+ S(g 3,o,2* 

f Q ; io 

9,o,o’ 

,40,30) + 

T I s(f 3 

,o,l’ 8 9,o,l ;l0 ’ 40 ’ 30) 

* S(g 3,o,l’ 

f 9,o,i ;l ° 

,40,30) ] | 

" C ll,2 

( S(f 5,o,o’ g 7,o,2 ; 


X 20,48,28) + S<g 5|Oi2 ,f 7jOjoi 20,48,28) ♦ £ [*(f 5>0>1 , 
x g ? ^20,48,28) * s <g 5 „ b , i“, 48,28) ] }, (A22b) 
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f U.o,2<°> 


f ll,o,2 ( °> 


8 11,o,2 (0) 


8 LI,o,2 (0) 


0 , 


f» o' + 0 

r ll,o,2 8 1 l,o, 2 


ai n ♦ •. 


( A22c) 
( A22d) 
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Table 2. Parameters of 

the elliptic cylinders 

for which results are given. 

a/b 

k 2 

E(k) 

5 

0.96 

1.05050 

10 

0.99 

1.01599 

25 

0.9984 

* 

1.00329 
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Table 3, The position in the boundary layer of the maximum of 

a 

a/b » 

5.0 

10.0 

25.0 

w/4 


2,2 

2.1 

2.0 

it 


1.9 

1.8 

1.7 

3* 


1.2 

1.1 

1.05 

6w 


0.9 

0.8 

0.75 


Table 4. 

Total 

phase shift, 

in degrees, of 

u i * 

1,0 

a 

a/b - 

5.0 

10.0 

25.0 

w/4 


11.2 

11.9 

12.0 

TJ 


30.0 

30.0 

30.0 

3W 


42.0 

40.0 

40.0 

6w 


42.0 

42.0 

40.0 
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Table, 5. Coefficient* in the power aerie* expansion of the potential 
velocity function, H(5)> for flow paat a parabolic 
cylinder. 


2n ♦ 1 


a 


2n 1 


l 


1/4 


3 


- l/2 2 .3! 


5 22/2 4 .5! 

7 -1168/2® .7! 


9 113536/2® .91 

11 -17521024/2 1 0 .ill 
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to 0. 74 with an interval of 0. 04. The body is an ellipse with 
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Figure 5. Contours of the second order steady flow, u 2 Q . Contours fro* 

-1.6 x 10~ 2 to 5.2 x 10“2 with an Interval of *4.0 x 1(T 3 . The 
labels are scaled by 10 4 . The body is an ellipse with a/b - 5 
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-1. 0 x 10” 3 to A. 8 x 10” 3 wlth an interval of 8. 0 x 10"A 
labels are scaled by 10 5 . The body is an ellipse with 
and o “ 6v. 







Ur - ■ 








gure 16. Contours of the second order steady flow, u 9 . Contours fros 

-7.7 x 10~ 3 to 4.2 x 10~ 3 with an Interval of 7.0 x lO” 1 *. The 
labels are scaled by 10 5 . The body is an ellipse with a/b - 25 
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Figure 17b. Contours of the phase, in degrees, of u_ the fundamental 

oscillation. Contours froa -11.2* to 0* with an Interval of 
0.7*. The body is an ellipse with a/b * 5, and a m »/4. 



ORIGINAL PAGE® 
OF POOR QUALITY 


e • *-«*-.* S« 



65 


oscillation. Contours fro« -28.0* to 2.0* with an Interval of 
2.0°. The body is an ellipse with a/b « 5, and o - 5. 
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Figure 19a. Contours of the amplitude of u n the fundamental oscilla- 

tlon. Contours from 0.0 to 0.96 with an Interval of 0.06. T1 
body is an ellipse with a/b « 5, and o - 3*. 
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Figure 21b* Contours of the phase. In degrees* of u_ « , the fund seen tel 

oscillation. Contours free -11.2* to 0.7* with an Interval of 
0. 7 *. The body Is an ellipse with a/b * 10, and a ■ t/4. 
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Figure 22a. Contours of the amplitude of u. the fundamental oecllla- 

tlon. Contours from 0.0 to 0.90 with an Interval of 0.05. T1 
body is an ellipse with a/b ■ 10, and o m a. 
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oscillation. Contours fraa -28. 0“ to 2.0* with an Interval of 
2. 0*. The body Is an ellipse with a/b ■ 10, and <j ■ t. 




i-1 I l ilt 1 1 1 1 1 1 1 1\| tl | Lll j 1 |/M » l 1 m i i i/i i I I I 


ft A 


ORIGINAL PAGE 13. 
OF POOR QUALITY 



•> 

£ 


d 

o 


S 

U 

8 

•a s 
8 

8 

■o 

§ 


o 

o 

3 

c 

0 • 
4J * 

a to 


J Sb 

a n 

. 3 g 

H S 
* * 
O O o 

<3 0\ rl 
o i 


•S 

3 

4J 

•H 


a > 

5 


S 

§ 

0 

O 

w 


O XI 

w *«5 

Sfi 

8 S 

M <U 
4-t 00 
3. 
to 

U H 
3 H 
O V 

4J 

o 0 
o n 
u 

00 


g & 

O T) 
•H O 
U A 


9 

CO 

es 

s 

Si 


f 

I 


74 



ORIGINAL PAGE IS 
OF POOR QUALITY 



U-4 

o 


* g 

u u • 

0 «l J5 

v u<n 

!*■ 
s r 

*« t> 

1 « 8 

3 v - 

• H 
I 

O 02 
< 3 O 

4J 0 

* e US 

MOW 

8 J© 1* 

m n 
Mi I 0 
0 a 

*° 8 -H 


0 2 a 

So" 


4 u id 
ft. a -h 
o 

0 o ►> 

3 *3 

• ^ 

*8 O 0 

g « £ 

h 0 

9 H 

O H . 

4J tH • 

a y o 
o o • 
U 0 CM 



75 



ORIGINAL PAGE » 
OF POOR QUALITY 



76 


Figure 24.. Contour, of the auplitude of S o l . the fund«ntal oaeill 

tion. Contour, fro 0.0 to 0.90 ’with an internal of 0.05. 
body is an ellipse with a/b - 10, and a - 6ir. 
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oscillation. Contours frosi -12. 0° to 0.0* with an interval of 
0.8*. The body is an ellipse with a/b ■ 25, and a ■ t/4. 
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harmonic of the free stream oscillation. Contours from 
-240.0° to 100.0° with an interval of 20.0°. The body is an 
ellipse with a/b = 5, and a = tt. 
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harmonic of the free stream oscillation. Contours from 
-260.0° to 80.0° with an interval of 20.0°. The body is 
ellipse with a/b - 5, and a - 3*. 
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Figure 33a. Contours of the amplitude of u 0 ^ 2 » the first harmonic of the 

free stream oscillation. Contours from 0.0 to 4.0 x 1G“ 2 with 
an interval of 2.0 x 10" 3 . The labels are scaled by 10 4 . The 
body is an ellipse with a/b * 10, and o * ir/4. 




harmonic of the free stream oscillation. Contours from 
-200. 0° to 140. 0° with an interval of 20. 0*. The body is 
ellipse with a/b * 10, and a * ir/4. 
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Figure 35a. Contours of the amplitude of u. «, the first harmonic of the 

O J 

free stream oscillation. Contours from 0. 0 to 7.2 x 10~ 3 with 
* an interval of 4.0 x 10 -1 *. The labels are scaled by 10 s . The 
body is an ellipse with a/b « 10, and cr - 3*. 
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Figure 36a. Contours of the amplitude of u the first harmonic of the 

c . °»2 

tree stream oscillation. Contours from 0.0 to 4.0 x 10~ 3 with 
an interval of 2.0 x 10~ 4 . The labels are scaled by 10 s . The 
body is an ellipse with a/b * 10, and a * 6*. 
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Figure 39a. Contours of the amplitude of u_ 9 , the first harmonic of the 

free stream oscillation. Contours from 0.0 to 6.4 x 10~ 3 with 
an interval of 4. O x 10“ 4 . The labetls are scaled by 10 5 . The 
body is an ellipse with a/b « 25, and a - 3*. 




stream oscillation. Contours from -260.0° to 60.0° with 
interval of 20. 0°. The body is an ellipse with a/b ” 25 
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Figure 40a. Contours of the amplitude of u 9 , the first harmonic of the 

°» z , 3 

free stream oscillation. Contours from 0. 0 to 3.4 x 10 3 with 
an interval of 2.0 x 10 -lf . The labels are scaled by 10 5 . The 
body is an ellipse with a/b “ 25, and a * 6s. 
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of -the free stream oscillation. Contours from -280.0" to 60.0 
with an interval of 20.0". The body is an ellipse with a/b * ; 
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Figure 41. Contours of the amplitude of u the fundamental oscilla- 

tion of the flow past the parabola with a * ir/16. Contours 
from 0.0 to 0.22 with an interval of 0.01. 












■abola with o “ 3*/4. Contours froa —3.0 x 10“ 3 to 2.4 x 
' 3 with an interval of 3.0 x 10 —l *. The labels are scaled by 
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of the flow past the parabola with a * ir/16. Contours 
0.0 to 0.22 with an interval of 0.01. 
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Figure 46b. Contours of the phase, in degrees, of u 0j ^> fundamental 

oscillation of the flew past the parabola with a = w/16. Contours 
from -8.5° to 0.0° with an interval of O.i*. 




tion of the flow past the parabola with a m x/4. Contours from 
0. 0 to 0. 23 with an interval of 0. 01. 
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oscillation of the flow past the parabola with a = t/4. Cbntours 
from -24.0° to 1.0° with an interval of 1.0*. 
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if the flow past the parabola w 
.0 to 0.22 with an interval of 





oscillation of the flow past the parabola with o ** 3w/4. Contours 
from— 34.0* to 2.0* with an interval of 2.0* 
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of the flow past the parabola with a “ 3ir/2. Contours 
0. 0 to 0. 22 with an interval of 0. 01. 
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Figure 50b. Contours of the phase, in degrees:, of u of the first harmonic 

. o, z » 

of the free stream oscillation of the flow past the parabola with 

a =r/16 3 contours from -260.0* to 140.0* with an interval of 

20 . 0 *. 
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Figure 51b. Contours of the phase, in degrees, of « o 2 * of Che firsC harraonic 

of the free stream oscillation of the flow past the parabola 
with a = ir/4, contours from to —240,. 0* to 100.0 with an interval 
of 20. 0 # . 
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ABSTRACT 


A solution of the Navier- Stokes equations is given for an incompressi- 
ble stagnation point flow whose magnitude oscillates in time about a con- 
stant, nonzero, value (an unsteady Hiemenz flow). Analytic approximations 
to the solution in the low and high frequency limits are given and compared 
to the results of numerical integrations. The application of these results 
to one aspect of the boundary layer receptivity problem is also discussed. 
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1, INTRODUCTION 

In this paper we give a solution to the Navier-Stokes equations for an 
incompressible stagnation point flow whose magnitude oscillates in time 
about a constant, nonzero, value. Apart from the intrinsic interest of this 
problem, its solution is the first step in the solution of one aspect of the 
boundary layer receptivity problem, that is, the determination of the magni- 
tude and form of the disturbance introduced into the boundary layer on a 
body by a perturbation in the free stream. The solution of this problem 
would permit the calculation of the initial amplitudes of the Tollmien- 
Schlichting eigenmodes and continuum eigenfunctions in the boundary layer 
and give a rational foundation to transition prediction methods. 

We have recently given the solution to the boundary layer receptivity 
problem within the context of incompressible, linear stability theory for a 
parallel shear flow (Salwen, Kelly, and Grosch, 1980; Grosch & Salwen, 1980; 
Salwen and Grosch, 1981). There is, however, one aspect of the boundary 
layer receptivity problem to which our parallel flow solution is clearly not 
applicable. If we consider the flow near the forward stagnation point of a 
body, a linearized parallel flow theory cannot be valid because the flow is 
intrinsically nonlinear and nonparallel. Although we may be able to use the 
boundary layer equations away from the stagnation point, the full Navier- 
Stokes equations must be used in the immediate vicinity of the stagnation 
point. Once a solution of the Navier-Stokes equations for the perturbed 
stagnation point flow has been found, it is possible to extend this solution 
away from the stagnation point and around the nose of the body. This is, in 
fact, what is normally done for the steady flow past an object. The Blasius 
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series (Schlichting, 1979, pp. 168ff) is the extension around the nose of a 
blunt body of the Heimenz soluton (Schlichting, 1979, pp, 95ff) to the 
Navier- Stokes equations near a stagnation point, In section 6, we discuss 
the analogous extension for our oscillating flow solution. 

We were stimulated to do this work by the papers of Glauert (1956) and 

Stuart (1966), There are a number of other relevant studies in this area 

including those of Lighthill (1954), Rott (1956), Lin (1956), Carrier and 

DiPrima (1957), Gibson (1957), Watson (1959), and Sarnia (1964), Riley 

« 

(1975) and Schlichting (1979, chapter 15) give comprehensive reviews of 
recent work in unsteady boundary layer theory. Lighthill (1978) has 
recently reviewed the current understanding of the phenomenon of acoustic 
streaming; i.e,, the generation of a steady flow by the Reynold’s stress 
due to an oscillating flow. 

2, EQUATIONS AND BOUNDARY CONDITIONS 

We consider the flow in the neighborhood of a stagnation point at 
(0,0) on a plane wall, with x the coordinate along the wall and y the 
coordinate perpendicular to the wall. We seek the solution (u(x,y,t), 
v(x,y,t)) of the two-dimensional Navier-Stokes equations which corresponds 
to the potential flow 


U * (U x/£)(l + e cos wt) 
0 

(1) 

V ■ - (U Q y/£)(l + e cos ut) 

(2) 
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in the Ear field. Here, U and V are the x- and y-components of velocity 
of the potential flow, U Q the velocity scale, % the length scale, e 

the dimensionless amplitude of the oscillation, and u> the frequency of 
oscillation. Defining a stream function 'Kxjyjt) by 


3tp 8if/ 

Ty* V cTx 


(2a, b) 


and substituting into the incompressible Navier- Stokes equations we find, as 
usual, that is the solution of 

( 5 T + 57 k- slsy 5 72 * ■ v72<72 * ) - < 3 > 

with v the kinematic viscosity and 


7 2 = iL + iL 

9 9 y^ 

We define dimers ion less variables 

K 3 x/A , n = (y/A)/R , t - u t/i, 

o 0 


(4) 


(5a,b,c) 


a Reynolds number, 


R * U J l/v 
o 0 


(6 a) 
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and a Strouhal number 


o . (6b) 

o 

The use of separate velocity and length scales, U q and Jt, while 
correct, is somewhat arbitrary in that they appear only in the combination 
U q /Z- in the potential flow, The scales which are intrinsic to this problem 
are the time scale of the base flow, T q * £/U q , the frequency, w, and 
the kinematic viscosity, v. We can define velocity and length scales in 
terms of T q and V, thereby setting the Reynolds number equal to unity, 

In any case, the dimensionless parameters which appear in the equations and 
boundary conditions are the oscillation amplitude, £, and the Strouhal 
number, r/ , and the results will be the same. We have chosen to use 
independent velocity and length scales, U q and £, because we intend to 
use the solution presented here as the basis for constructing a solution to 
the problem of an oscillating flow past a blunt object with length scale l. 

If we set, in analogy with the Hiemenz solution (Scnlichting, 1979, pp. 
95ff), 


- (0 A//R ) 5 F(n ,t) , 


then equation (3) reduces to 


(LL) 2 - 

3q 3x3n 3n 


3 2 F _ 3 3 F 
dvF Sq " 3 


•] ■ 0 , 


(7) 


( 8 ) 


HO 



which may be integrated to yield 
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5ii_ + <ii* 

3x3n 3n 


3 2 F _ 9 3 F 

3n 2 3n 3 


G(£ ,x). 


( 9 ) 


Because the left-hand side of equation (9) is independent of $, it is 
clear that the "constant of integration", G, only depends on the 
dimensionless time, t, in order that the x and t dependence of the 
pressure in the far field agree with that of the potential flow, we must 
have 


(U A ) 2 G + (U /A ) (~ew sin wx) + (U /l) 2 

X 3q x 3 x 0 0 


x (1 + e cos w t) 2 . 


Therefore, F(n,x) is the solution of 

3^_F_ + ^3_F^ 2 _ F 3 2 F _ 33 f » i 

3x3n 3n 3n 2 3n 3 

+ y e 2 (1 + cos 2 ot ) 


e(2 cos ox - o sin ax) 


(10) 


( 11 ) 


with the boundary conditions 


F(0,x) * 0, 


(12a) 


On-o * °» 


(12b) 
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«srr + 1 + cos or as n + ”, (12e) 

an 

3, METHOD OF SOLUTION 

To solve equation (11), we will expand F(n,t) in a Fourier aeries in 
Che dimensionless time, T. The coefficients of the expansion are functions 
of n and the parameter e and each of these will be expanded in a power 
series in e, It is easy to see, from the form of (11), that the x~ 
independent term in the Fourier series contains only even powers of e. 
Therefore, we look for a solution of the form 

F(n,r) = l U 2m f <n) * J e 2m+k [f <n) =° 3 tot 

m*o * k*l * 


+ *2mA ( ' l) Si " k0TlK 


(13) 


To find the equations obeyed by the f 2m k and g 2m ^ we substitute 
(13) into (11), collect the coefficients of like terms in the Fourier 
series, and set the coefficients of the successive powers of e equal to 
zero. We find, first, that f is the solution of 


d 3 f d 2 f df 2 

+ f - ( — 2 jl°) +i = o (14) 

dn 3 0,0 dn 2 dn 


with 
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(15a,b,c) 


Wvpi ffj 

w prru o v-v 


£ (0) " V to) * 0, f +■ 1 m n + «, 

o,o x ' fc o,o * o,o 1 


As was expected , f , the T- independent, C- independent term in the 

w | w 

series, is the Hieraenz solution for the steady stagnation point flow. 
Next we define the operator L by 


L s j£- + f. 

dn* 


d2 - 2f 1 1- + f" , 

dn* °»° dn ” 


(16) 


with primes denoting differentiation with respect to h. Then it can be 

shown that the equations and boundary conditions for the {f„ , } and 

zm, <c 

‘ Vk 1 

are, for k 2. 1, 


Lf„ , - ko g 1 0 , ■ P. , , 

2m, k 2m, k 2m, k 


Lg. , + ko f ' * Q rt , 

2m, k 2m, k x 2m,k 


(17b) 


J 2n,,k <0) ’ 8 2m,k (0) * £ '2m,k (0) ’ 8 '2 m> k <0) " °> 


f' + 6 6. . , g' . * 0 as n + °°, 

2m, k m,o k,l 2m, k 


(17c) 

(17d) 


and, for m > o and k ■ 0, 


L ^2m,o ^2m,o* 


(18a) 
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H. 


..ikinl DAft E IS 


£ 2 m ,o<°> ’ £ Vo<°> ■ »■ 


f' + 0 as n + 
2m, o 


(18b) 
< 18c) 


The (P„ . ) , (Q 0 . } , are linear combinations of products of the 

Zin | k zm,K 

( £ } , { g 0 ) and their derivatives, where a < k and 0 S. r £ ro or 

2r,s * °2r,s 

s ■ k and r < m, and are given in Appendix A. Therefore, these equations 
can be solved sequentially. The equations can be integrated numerically 
quite easily, although care must be taken to control roundoff errors. The 
results of these integrations are given and discussed in Section 5. 

4. LOW AND HIGH FREQUENCY APPROXIMATIONS 
4. 1 The Low Frequency Limit 

As o+0, it is expected that the solution will approach a quasi- 
steady solution. It is straightforward to show that the quasi- steady 
solution, 


F(n ,t ) * (1 + e cos or ) 1/2 f ([ 1 + e cos ar] I/2 n), (19) 

O f 0 

satisfies equation (11) and the boundary conditions (12) to 0(o), 

It is also easy to show that this quasi-steady solution is consistent 
with the expansion given in equation (13). If F, as given in (19), is 
expanded in a Taylor series in e , we find that 


f - l i£-52L2li-s (n ) 

n*o 2 n n! n 
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*-’3 ft . 

« I' / 


00 “ 2ntf-k S„ , , (n) 

1 <2 - «„„) I — - - ■ CO. to»T, 


k*o 


ko / t '4 
m*o 


ra! (m + k) ! 


(20) 


with 


s (n ) * f (n ) , n “ o 

n 0,0 


n- 1 n- r- 1 


n-f^Oi) + l (-1) (2n- 5 -2)1 

o,o „n- r- 1 ,/ v, n o.o 

* r^o 2 r!(n-r)I 2 


n > 1, 


( 21 ) 


and 


[<*) 

o.o 


d r f 


0,0 


(2 2 ) 


dr ( 1 


Equating the coefficients of e in (21) to the corresponding 
coefficients of e in (13) we find that, to order a, 


‘2m, 0 2 4m (ra!) 2 S2ra, 


(2 3a) 


£2m ' k ° 2^ m+2k “ ^(m+k) Im! S2m¥k ' 


(23b) 


2 = 0 

°2m, k 


(23c) 
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Taking Che a + 0 UmiC in the differential equations (17) and (18), for 

the { f„ .} and (g„ and substituting the low frequency approxima- 

2m ) k 2m ) k 

tions (23), it can be shown by induction, after some lengthy but straight- 
forward calculations, that these approximations satisfy the dif ferentional 
equations to 0(a). Therefore, the quasi-steady solution is the same as the 
expansion in the low frequency limit. 

4.2 The High Frequency Limit 

In the high frequency limit, a -*■ ", it would be natural to look for 

approximate solutions for the { f n .} and {g } in the form of a power 
rr 2m, k 2m, k 

series expansion in a"** , However, it is clear from the form of the 

differential equations, that this expansion would be non-uniform because the 

highest derivatives would be multiplied by the small parameter a" 1 and 

thus would vanish as a ->• 00 . 

We therefore rescale the equations, defining a new independent 
variable, an inner variable, 

z => (a/ 2) 1/2 n * (w/2V) 1/2 y, (24) 


with the length scale (2v/w) 1/2 , that of Stokes's second problem 
(Schlichting, 1979; pg . 75 and Chapter 11). In this limit we assume that 
there is an inner expansion of the form, 


o.o 


(n) + 


l 

nt= 1 


2m 


f 2m,o (z) 


00 00 

+ l l e2ro+k [£ 2m C0S kaT + ^2ra k^ sin k0T ^ 
m«o k=l 
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It is easy to see that the differential equations for the (f 2m 
and {g„ . } are just equations (17) and (18) with a transformation of 

fcin | K 

variables from ri to z. It can be inferred from the form of the equations 
that f„ , and g 0 . are 0 ( ( 2 / ct) 2m+k '' 1 / 2 ) i n t h e high frequency limit. 
We then solve these equations, retaining only the highest order terms in an 
expansion in powers of ( 2 /o). 

We find 

f , (z) * (2/a) /2 [z + ~ e~ z (cos z - sin z) - - 5 -], (26a) 

0,1 z z 

g© x (z) ■ ( 2 /o) [■— e z (cos z + sin z) -y]. (26b) 

a Stokes shear layer flow caused by the 0(e) part of the far flow field. 

In the high frequency limit it is decoupled from the 0(e°) steady outer 

flow. This is, of course, a familiar result in time dependent boundary 

layer theory and has been derived and discussed by Carrier and DiPrima 

(1957), Stuart (1966), and Riley (1975), among others. 

If we let 6 and 6 , be the boundary layer thickness of the 

0,0 o,l 

steady flow and the Stokes shear layer flow, respectively, then, for 

a » 1 , 


6 ,/6 ** 2 (2/o ) 1/2 « 1. (27) 

o,l 0,0 

Next, we solve the equations for the 0(b 2 ) oscillatory flow and find 

that 
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f „(z) ■ (2/0 ) 3/<2 {^Z[e >/r ^' Z (sin fl z + cos fl z) - l] 
0>l 16 


1 -z . . 

+ ^. ze sin z> 


(28a) 


g - (z) = (2 /a ) 3 /2 {— [ e ^ Z (sin /I z + cos fl z) - l] 


3 o,2 


16 


- i ze cos z} (2 8b) 

Again we have a Stokes shear layer, decoupled from the steady flow and, for 
large a, confined to a thin layer imbedded within the steady boundary 
layer . 

Proceeding next to the 0(e 2 ) steady streaming flow component, we find 
that the most general solution which satisfies the boundary conditions at 
z = o, is 

f„ = (2/a ) 3 1/2 (-L-) [ 13 - 6z - e - 4e Z (3 cos z + 2 sin z) 

2 ’° 16 


- 4ze Z sin 


z + Dz 


(29) 


where D is an arbitrary constant. This secondary steady streaming flow is 
identical to that found by Stuart (1966) using the boundary layer equations. 
Stuart's small parameter, a, is the reciprocal of o and Stuart's expan- 
sion is in powers of a, while ours is in powers of (2 /a). Therefore, 
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p AGB is 
QUALITY 

f , with D * 0, is equal to 1/2 of the x given by Stuart (see Stuart, 

2,o s 

1966; eq. 2.1L), Our solution and Stuart's satisfy the boundary conditions 
at the wall but do not satisfy the outer boundary condition because 

df 9 n 

(2 /a ) + o Dn as n + » (30) 


ORIGINAL 
OF POOR 


and does not vanish, as required by the far field boundary condition, for 

any value of D. This is explained by Che fact (pointed out by Stuart) 

that there is an outer boundary layer, thick compared to the Stokes layer 

thickness, 6 ,, but thin compared to the scale of the body, within which 

o,l 


this secondary steady streaming flow decays to zero. 

Since the o + 00 limit is non-uniform and z is an inner variable 

A A 

all of the f 1 s and g' s are inner solutions and could be expected to 


UULC1 OU JLLl L XUUO • 


require matching to appropriate 

g , f , and g , however, are also outer solutions and, unlike 
o, 1 o,2 o,2 

f , do not require matching. 

2,o 


> 


Outside of the Stokes layer, the interaction between the secondary 

A A t 

streaming flow, V , and the oscillatory potential flow, f ' 1 , is 
unimportant. In the case studied by Stuart, the dominant non-linear 

✓N 

interaction is that between f' and itself. Since the velocity of this 

2 ,o 

flow at the top of the Stokes layer is 0(U Q /a), the corresponding Reynolds 

number is R = (U /a)(£/v) = R /a. From this, Stuart concluded that the 
so o 

thickness of the outer layer is 0(^R“^^) = 0(&R -^2 CT 1/ 2) ^ w hich is 

so 

larger than that of Stokes layer by a factor of the order of a. 
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For the case under consideration here, the situation is completely 

changed by the existence of the large steady flow, f^ Q with boundary 

laver thickness 6 ■ 0 (cr 1/2 6 The dominant non-linear interact! 

J 0,0 0,1 

of f 0 in the outer layer is with f* and, consequently, it falls off in a 
2,o o,o 

distance of the order of 6 (instead of o 1/2 5 ). The appropriate 

0,0 0,0 

outer variable is therefore n =* y/ 6 


ion 


o.o 


If we express f . and g . in terms of n, substitute them into 
r o , i o,l 

the differential equation for f^ , and let a + °°, it can be seen that the 

A 

outer solution, f 9 , is the solution of 

, 0 


L (f ) a (i + f 2f 1 — + f" ) f„ - 0, 

v 2,0 , 0 o.o , 0 o.o ... 0.0 2,0 

1 dn u dn- ' dn 


(31) 


with the outer boundary condition 


f ' * o as n 08 . 

z, o 


(32) 


Using the asymptotic expansion for f , 

0,0 


f ~ n + A, 

. 0,0 


(33) 


with A a constant, it can be shown that the general solution of (31) is 


f„ = c f' (n) + c,h , ( n) + c 7 h 9 (n) 

2,0 o o,o 1 1 ' L L 


(34) 


where, as n », 


2 / 

h' ~ (1 + 0 2 ) / e" Y ,2 dy - 0e 

1 S 


-3 2 /2 


(35) 
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0 + y 0 3 , 


(36) 


with 


0 * n + A. (37) 

The outer boundary condition requires that C2 be set equal to zero. 
This leaves three arbitrary constants, C q , , and D, in the inner and 

outer solutions so we can match the inner and outer solutions and their 
first two derivatives. Matching the inner limit of the outer solution and 
the outer limit of the inner solution shows that 

C Q « OCa- 1 ), C x - 0(o- 3/2 ), and D - Q(o~ 2 ). (38a, b,c) 

i y a 

The inner solution was obtained by expanding in powers of o' and 
retaining only the lowest order terms. Consistent use of this approximation 
requires that we set C. and D to zero and match the first derivative. 

We find that 

C --3/ (ABO). (39) 

o 

where B * f M (0). The composite solution for f is 

0,0 2,0 

?2 0 * (2/o) 3/2 (-^g.) [ 13 - e 2Z - 4e Z (3 cos z + sin z) 
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- 4ze 2 sin z] — 6' (n), 

4Ba °*° 


(40) 


which satisfies the boundary conditions at z equal to zero as well as 

z ln . The thickness of the inner layer of the steady streaming flow is 

^ , while that of the outer layer is q which is much larger than 

5 . for o *► <*>. 

o,l 

We have also calculated the high frequency approximations for the 

0 (e 3 ) functions, for f Q g, g Q 3 * ^ 1 and g 2 1 * Tltese are gi ven i° 

Appendix B, The components of the amplitude of the 0(e 3 ) portion of the 

second harmonic f 0 and g 0 are driven by the interaction of the 0 (e) 
fundamental (f , ,g . )and the 0(e z ) first harmonic (f _,g _). The 

0 | 1 0)1 0 , Z 0 5 4 

inner expansions for f „ and g -, given in Appendix B, are also outer 

o j j OjO 

expansions • 

This is not true for f 2 ^ and g 2 ^ , the components of the amplitude 

of the 0(e 3 ) part of the fundamental. They are excited by the interaction 

of (f ., g ,) with (f 0 ) as well as (f . ,g^ .) with (f „,g «). 
o,l o,l 2 ,o o,l' 0,1 o ,2 0,2 

The inner expansion of the in-phase component, f 2 satisfies the outer 

boundary condition, but inner expansion for the out-of-phase component, 

g , does not. The outer expansion for n is the same as that for 

2,1 ’ 

f . The composite expansion for g„ . is found in the same way as that 

1 0 £ y L 

for f 9 and is given in Appendix B. 

4 j O 

5 , RESULTS 


We have numerically integrated the differential equations for the f's 

and g' s with 2m + k < 3 over a wide range of value?/ of a, Some typical 

/ 
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results at a moderately low frequency, a » ir/4, are shown in Figures 
1-4, where the solid curves show the results of the numerical integrations 
and the dashed curves are obtained from low frequency approximations for the 
f's given in section 4.1. The g's are identically zero in the low fre- 
quency approximation. 

It can be seen from the results shown in these figures that at low 
frequencies there is no Stokes layer; the boundary .layer thickness of the 
various components are generally equal to or greater than that of the steady 
flow component, f . The steady streaming component, f' , is quite 


small compared to the mean flow, even for e - 1. 

At ct ■ o the g's, the out-of-phase components, are identically 
zero. At small, but non— zero, frequencies the low order (m s) in-phase 
components, the f's are larger than the g’s. As the order increases, 
however, the magnitudes of the f's and g's tend to equalize and 
decrease. 

It is also apparent from an examination of Figures 1-4 that the low 

frequency approximations are reasonably accurate even at 0 * tt/ 4. In fact, 

the difference between the numerical solution and the approximation for 

f is so small that it is not apparent in Figure 2. Taking into account 

o, 1 


the fact that the higher order terms, which have the largest deviations from 
the low frequency approximations, have very small magnitudes, it is clear 
that the low frequency approximations, equation (19), is reasonably accu- 
rate, even for Strouhal numbers as large as r/4. 

Figures 5-9 show the f's and g's for 2m + k < 3 at a moderately 
high Strouhal number, o «* &i , as obtained from the numerical integrations 
(solid curves) and the asymptotic approximations (dashed curves) . The top 
of the Stokes layer is at z * 4.6 and this corresponds to n " 6.5 
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Therefore, the top of the Stokes layer is at n * 1,3 for a " 8tu The 

secondary streaming flow (see Fig, 5) extends far beyond the Stokes layer, 

The variation of the fundamental component (see Fig, 6) is confined to the 

Stokes layer. The amplitude functions of the higher harmonics tend, like 

f* , to extend beyond the Stokes layer because they are also driven by 
Reynolds stresses due to f Q Q . 

The secondary steady streaming flow, f'„ , is considerably smaller 

1 1 o 

at high frequencies than at low frequencies. In addition, the nature of the 
secondary flow changes as a increases; at a small o, the net secondary 
flow is positive while, for large a, it is negative. As discussed by 

Stuart (1963), it is known from experiments that this effect occurs in the 
steady streaming flow generated by a circular cylinder oscillating along a 
diameter, 

8F 

The tangential velocity component, u, is proportional to that 

is to the { f ^ g^ . in Figure 10 we present some of the results of 

the numerical integrations; a plot of the variation of max | f^ jJ and 

max I gi . I as functions of a. The maximum of f' is 1.0 for all o 

|0 2m,k' o,o 

and the absolute maximum of f . is about 1.069 at a " 17 ; the maxima of 

o,l 

f and f , are not plotted in Figure 10. We can conclude, from the 
o,o o,l r u 

results shown in Figure 10, that the high frequency estimate, {| f^ k l, 

1 82m, k I ^ * °^ a 2m k+1 ^> is quite good. 

6. DISCUSSION 

The applied far-field flow consist of a time-independent mean flow and 
a fundamental with a cosmt time-variation, Nm-linear interactions result 


154 


in the generation of components at ail multiples of the fundamental fre- 
quency and modifications of the mean flow and fundamental components, 
though, at moderate e(£ 1), the fundamental and first harmonic tend to 
dominate. 

At low frequencies, our results are well represented by a quasi- steady 
approximation (19), which has the same form as the steady stagnation- point 
flow except that the amplitude and scale vary with time, Somewhat surpris- 
ingly, this approximation is quite good for a dimensionless frequency, o, 
as large os ir/4. 

♦ 

At high frequency, the viscous boundary layer corresponding to the 

oscillating components is largely, but not completely, confined to a Stokes 

layer of thickness (2/o)^^ times that of the steady boundary layer. For 

these large values of o, the inner asymptotic approximations are solutions 

of differential equations which are independent of f 1 and are, there- 

0,0 


fore, decoupled from the mean flow. For a number of components, the inner 
expansions are also the correct outer expansions; these components are, 
therefore, totally decoupled from the mean flow to lowest order in 0 and 
are, in fact, identical to the corresponding solutions found by Stuart 
(1966) for a purely oscillatory flow past a body. 


Not all of the components are decoupled from the mean flow in the high 

frequency limit. The high frequency expansion is a non-uniform asymptotic 

expansion and thus, for certain components, in particular the steady second 

order flow, f' , an outer expansion, matched to the inner expansion, is 

2,o 


needed for a uniformly valid approximation. The results given in Section 
4.2 show that such composite expansions can be found for 0 >> 1; they 
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satisfy, to lowest order in er, the differential equations and all of the 
boundary conditions. We believe, although we have not carried out the cal- 
culations, that this is also true for an expansion to any order in o. The 
results of the numerical integrations, shown in Figures 5-9 are consistent 
with this belief. We conclude that, both for this flow and for the oscilla- 
tory flow studied by Stuart (1966), there are second, and higher, order 
steady streaming flow components that do not extend to infinity but are 
confined to a layer adjacent to the solid boundary. This layer is much 
thicker than the Stokes layer but much thinner than the length scale of the 
body. 


There are other aspects of Stuart's solution which can be compared w: 

II- - J i. t- — . — — . fu d AM J D « m 1% m ** Aw f.«U 1* rtK 

ulira» nc Uocu uiituc pac oiiicuci a > v* t hj uiiu *\ iu mao oyiai.^vuj w»*vn 

8 

are related to our Strouhal number, a, and Reynolds number, R q , by: 


a « U o /(tuJZ, ) * a- 1 , 

S + v/(w£ 2 ) « 1 /(or q ), 

R ■ U 2 /(wv) - R /a . 
so o 


Stuart used the boundary layer equations and an expansion for a + o; we 

have used the Navier-Stokes equations and, in addition to the numerical 

solutions, an expansion for the high frequency limit, o + It is clear 

that Stuart's expansion, and ours, are high frequency approximations. 

Stuart assumed that 6 was small and R large in order to justify the use 

s 

of the boundary layer equations. In the high frequency limit a + 00 , and, 
with R fixed, a, 0, and R are all small. We do not need to make 

O Jit s 

any assumptions concerning the magnitude of R q or R g because we used the 
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Navier-Stokes equations. We can determine the f 'a and g’s by numeri- 
cally integrating the ordinary differential equations; however, the quasi- 
steady solution for low frequencies and the high frequency expansion are 
useful analytic approximations. 

We have not attempted to determine the radius of convergence, in e, 

of the series for F(n,T,e), equation (13). However, some observations can 

«*2m-k+l /2 

be made. We have found that f 2m k and g 2m k are 0(o ) for 

large 0 . Therefore, max { | f^ k i , |g' 2m the ain P litucle functions for 

u, are O( 0 **^ m, *^ + ^) , Thus the series for u is in terms of (e/o)^ m+ ^, 

for large 0, and this suggests that it converges for e/o < 1 and 

converges rapidly for e < l, On the basis of the results presented in 

figure (10) we conclude that the high frequency bound on max {|f' 2ro k l » 

|g' |) is valid at all frequencies and, therefore, that the series 

converges rapidly for e < 1. Of course the convergence will be slower at 

low frequencies but we have shown that, for 0 ~ ir/4, the quasi-steady 

solution is an accurate approximation. 

In the introduction we suggested that the solution to the problem of 

oscillating stagnation point flow was the first step in the solution of one 

•» 

aspect of the boundary layer receptivity problem. A few years ago, Morkovin 
(1978) reviewed the rather rudimentary state of knowledge of the dynamics of 
boundary layer receptivity. He identified four general classes of free 
stream disturbances which might generate Toilmien-Schlichting waves in the 
boundary layer. These are; vorticity fluctuations, sound, entropy 
disturbances, and unsteadiness in the mean flow. In the incompressible 
limit, there can be only vorticity fluctuations and unsteadiness. Morkovin 
argues, however, that a temporally oscillating incompressible flow is a 
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reasonably accurate approximation for a sound wave impinging normally on the 
nose of a body if the wavelength of the sound wave is much greater than the 
radius of curvature of the nose, 

The solution of the Navier-Stokes equations given here is the solution, 
in the neighborhood of the stagnation point, of the receptivity problem for 
a simple unsteady flow. In order to interpret this solution in terms of a 
stability model it is necessary to extend it around the nose of the body 
into a region where the flow is, at least, quasi- parallel. If we can assume 
that the tangential component of the potential flow at the edge of the 
boundary is of the form U ■ U q H( 5) (1 + e cosat), with H(£) having a power 

series expansion in 5, the distance along the bod- from the stagnation 
point, we can use an unsteady variant of the Blasius series to solve this 
problem. For the firBt term we use the full two-dimensional Navier-Stokes 
equations and the solution is that given here. For the subsequent terms, we 
use the boundary layer equations. We have carried out this calculation 
procedure up to and including the terms in S 11 and we are now applying it 
to bodies which are elliptic or parabolic cylinders. We hope to report the 
results of these calculations at a later date. 

This work was supported, in part, by grant NAG- 1-96 from the National 
Aeronautics and Space Administration. 
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Appendix A 


The expressions for P 2m Q 2m and R 2m q in (17a), (17b), and (18a) 
are simplified by use of the notation 

M( f ,g) ■ f"g - 2f'g' + fg". (Al) 


In terms of M, 


P ok ' -X f8S U + 26 k2 * (1 - S kl> X [ M( t 0jr .f 0>k - r > 


- M(g ,g . «. ]} for k £ 1, 

°o,r ’°o,k-r) J ’ 


(A2) 


ra -1 


^2m,k 4 | g ^ Q ^^2(m-s) ,o ’ ^28 ,o 


m k-1 

• „> * <! ' V l I 


3=0 r=l 


X ^ M ^2s ,r ’ ^2(m-s) ,k-r^ M ^ s 2s ,r ’®2(ra-s) , k-r^ 


m-1 m-l-s 

+ 2 ^ ^ f M(f 2s,k+ra-s-s' jf 2s' ,m-s-s ' ^ + M ^ s 

3=0 3 ~0 


2s , k+m-s-s 1 ’ 


®2s ,m-s~s ' ^ | 


for m > 1, k > 1, 


(A3) 


k-1 


Q ok = 2 } 2a6 kl 


<X “ 6 kl> M(f o,r»*o,k-r> f0r k ~ 1 ’ (M) 
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‘2m, k 


m ~l v m k-1 

\l M[f 2U- s ,.S 2Slk ) * (l - V .I. J, 


8*0 r“l 


m-1 m-l-s 

X ^ f 2s,r ,8 2(m-s) ,k-r^ + j ^ f 2s, m-s-s ' 


s*o s'*o 


8 2s 1 ,k+m-s-s' * ^^2s 1 , k+m-s-s 1 * 8 2s ' , m-s-s ' ^ 


for m-1, k-1, 


(A5) 


and 


m-1 


R 2m,o " -4 Z4 „,l + - S „,> l 


ml ' 


2s, o’ 2(m-s),o 


) 


m-1 m-l-s 


a q'=o ^ M ^2s , m-s-s 11 f 2s ', m-s-s ' ^ + M ^ 8 2s , m-s-s ' , 

v v 5 *"U 


* S 2s' , m-s-s ' ^ 


for m-1. 


( A6 ) 


Here, is the Kronecker 6 symbol. 
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Appendix B 

High Frequency Approximations to the CXe" 1 ) Functions 
The approximations to f Q ^ and 8 0 3 f° r 0 * “ are; 

5/2 

f (a) - (2 A?) {[9 - 8/1 + 3 (2 + l)) + e~ Z [(/7 

0 , J 

2 2 
-6 + 4z - 4z ) cos z + (6 - /T - 4z ) sin z] 

+ e ^ Z (4 + 2/T - 8z) cos /2z - ( 4 + 2 / 2 ) sin /2z] 

- / 3(/"2 + l) e ^ Z ( cos /3z - sin /3z) + ( 5/2 - 7) 

x e^"^" + ^ Z [ cos (/T +l) z - sin ( 2 + 1) z] } /128, (b 1) 

5/2 - z r / - 

g (z) * (2/a) {[9 - 8/2 + / 3 (/~2 + 1 )] + e [(/2 

0)3 

2 2 
-6 4z ) cos z + (/T - 6 + 4z - 4z ) sin z] 

+ e ^ Z [(4 + 2 / 2 ) cos flz + (4 + 2/T- 8z) sin /2z] 
-/I 

- + l) e Z ( cos /*3z + sin /3z) + ( 5/T - 7 ) 
x e - (/Tt 1)z [ CO s (/T + l)z + sin [Si * l)z]} /128. 


(b 2 ) 
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, # A t 

The approximation to . for o *► » is; 

l y 1 

5/2 

f 2 : ( z) - (2/a) {(349- 335/7)/ 19 20 + e [( 60/7-213) /9 60 

+ ( (5z/64) (z-l)/64) cos z + ((93 + 15/7) /240-lOz 
x (z + U/64) sin z] + e ^ Z [(7/7- 2 + 4z) cos flz 
- (7/7" + 2) sin /7z] /64 + (13/160)e"" 2z + e "^ +) - )z 
[(/7-5) cos(/7- l)z - (/7 + 5) sin (/Y- l)z]/384 
+ e ^ Z cos z/3 20} . 


The composite approximation to 



for a +■ 00 is: 


5/2 

g 2 x ( z ) = (2/a) ( - ( 1513 + 125/7) / 19 20 + e' Z [(558- 45/7] /9 60 

+ ( (3z (6-z)/64] cos z + ( (177 + 45/7) /960 + 
x z(3 L + 9z)/64) sin z] + e ( 7 / 2 " + 2) cos v7z 

— 9 7 

+ (7/7 - 2 + 4z) sin /2z]/64 + e [(13/80 + z/16)] 


(b 3 ) 
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with 


+ e ^ + ^ z [(5/2) cos(/2 - l)z - (5 - /2) 

A 

x sin (/*2 - l)z]/384 + e Z sin z/320} 

2 

+ ( 2 /a) 3f 1 (n)/(8B), 

0,0 

B - f" (0) . 

0,0 


(b4) 


(b5) 
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Figure 2. The 0(e), f^ and and 0(e 3 ), f£ ^ and g£ 

amplitudes of the fundamental component of the oscillating flow 
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the results of numerical integrations and the dashed 
the high frequency approximations. 







